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INTEGRATZCM OF A CODE FOR AEROELASTIC DESIGN OF 


CONVENTIONAL AND COMPOSITE WINGS INTO ACSYNT, 

AN AIRCRAFT SYNTHESIS PROGRAM 

By Joseph Mullen, Jr. 

Nielsen Engineering & Research, Inc. 

SUMMARY 

The integration, correlation, and documentation of the program for 
Wing Aeroelastic Design (WADES) of conventional and composite wing struc- 
tures for use with the aircraft synthesis program ACSYNT is described. 

A comparison of program estimates of wing weight, material distribution, 
structural loads and elastic deformations with actual Northrop F-5A/B data 
is presented. Correlation coefficients obtained using data from a number 
of existing aircraft are computed for use in vehicle synthesis to estimate 
wing weights. 

The modifications necessary to adapt the WADES code for use in the 
AC.sYNT program are described. Basic program flow and overlay structure 
is outlined. An example of the convergence of the procedure in estimating 
wing weights during the synthesis of a vehicle to satisfy F-5 mission 
requirements is given. A description of inputs required for use of the 
WADES program is included. Possible extensions and modifications of the 
structural model and analysis methods are identified where improvements 
in overall weight prediction and correlation with existing aircraft may 
be obtained. 


INTRODUCTION 

From 1972 to 1974, under the sponsorship of the National Research 
Council, structural optimization techniques were developed for the design 
of simplified conventional and multilayered composite wings for strength, 
stiffness, frequency, and flutter requirements. A computer program for 
wing aeroelastic design (WADES) was generated as the result of that 
investigation. The desirability of incorporating this capability into 
aircraft synthesis was identified so that the full impact of advanced 
concepts could be studied. Under Contract No. NAS2-0558* to NASA/Ames 

Research Center, Nielsen Engineering & Research, Inc. (NEAR) was funded 
— 

Technical Monitor; Dr. G. N. Vanderplaats . 






to incorporate the WADES program into ARC a ACSYNT program for vehicle 
synthesis. This is a final report summarizing that work. 

The primary purpose of this report is to present the results of the 
correlations and of the integration of the WADES program as a module in 
the ACSYNT program. A user's guide for the operation of the program is 
also included. Detailed comparisons of estimated weights, material distri- 
butions, and general assumptions with those of the F-5A/B wing are contained 
in this report. Correlation of the progrsr:.-e8timated weights with a 
broader group of U. S. fighter aircraft wing weights is also summarized. 

The program integration into ACS/NT was completed, and the results of a 
sample vehicle study are shown. 


SYMBOLS 


[A] 

{BI 

b/2 

^ig 

g 

(g) 

GLR 

®xy’®x2* ®yz 
h 

[Kr) 


matrix containing the aerodynamic influence coefficients 

matrix containing the unsteady aerodynamic influence coeffi- 
cients as derived for piston theory 

semispan 

polynomial function describing the shape of the wing camber 
surface 

polynomial function describing the wing semi -depth distribution 

arthotropic moduli of elasticity in x- and y-directions, 
respectively 

force on the landing gear on ground impact 

acceleration due to gravity 
vector of design constraints 
gross lift required 

orthotropic shear moduli of elasticity in respective coordinate 
plane 

altitude 

reduced stiffness matrix 
free-stream Mach number 
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;] conslst«nt mass matrix 

My banding and torsional moments about structural axis 

M* banding and torsional moments along airplane x,y reference 

^ axes 

load factor 

distributed aerodynamic pressure loading 

fraction of lift on wlngi at landing as fraction of total 
weight of vehicle 

free-stream static pressure 


distributed Inertial loadings due to skin, core, and fuel 

fraction of volume of structural planform available for fuel 

work equivalent load vectors due to aerodynamic, distributed 
loadings, and concentrated mass leadings 

stress gross weight 

wing planform area 


»n) 


f 


in 


»n) 


dy 


polynomial function describing the thickness of the skin 
over the wing planform 

wing depth to chord ratio 

"effective” skin thickness including distributed thickness 
of spar caps 

skin thickness of wing cover sheets only 

displacement vector 

shear load normal to wing surface 

function describing the transverse deformed shape of the 
wing 

total weight 

total weight, all components located in the body 


weight of concentrated mass 


el 


weight of fuel located in wings 


WOTO 

Wi 

ig 

'"^wing 

x,y,2 

x/c 


X V 
cm* -'em 


eg 


a 


A 


0_ ts« 

LE* TE 


G , a 

X* y 


*^xy * ^xz * *^y z 


U) 


vehlcl* gross weight at takeoff 
weight of main landing gear 

weight of wing 

spatial coordinates of wing 
local streamwlse fraction of chord 
coordinates of concentrated masses 

location of total weight as a percent of mean aerodynamic 
chord 

function describing the Internal structural rotations about 
the y-axls 

root angle of attack of wing 

function describing the Internal structural rotations etbout 
the x-axls 

angle of sweep of wing quarter chord 

nond'imenslonallzed coordinates of the wing; ^ » y/R and 
T • y/SPAN 

leading-edge zmd tralllng-edge sweep anglus for wing 
In-plane stress components In x- and y-dlrectlons 
shear stress components In respective coordinate plane 
dynamic frequency 
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PROGRAM METHODOLOGY 

The program used in the following study is the computer code for Wing 
Aeroelastic Design (WADES) developed under the sponsorship of the National 
Research Council. The program was developed for the preliminary design of 
conventional and multilayered composite aircraft wings to satisfy strength, 
stiffness, dynamic and flutter requirements. It models the structure of 
the wing as an equivalent orthotropic plate. The skin material distribution 
airfoil depth, and internal structure (core) are approximated by polynomial 
functions thac are continuous over the planform of the structure. The 
various static, dynamic and flutter analysis are performed using a Ritz type 
analysis with assumed polynomial modes. The program contains subsonic and 
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•upvrsonlc aerodynamic and Inertial static loads. Supersonic piston theory 
unsteady aerodynamics are used for flutter calculations. 

The design algorithms used in the WADES program are based on the mathe- 
matical programming technique known as the Method of Feasible Directions. 

The program employs either a direct search algorithm for the combined strength 
and dynamically constrained wing design, or a more efficient iterative proce- 
dure which uses a sequence of nonlinear approximate designs to converge on 
the least-weight design for strength constraints only. In each of these 
search techniques the parametric nature of the constraints on strength and 
minimum gage is converted to explicit form by evaluation only at discrete 
points on the structural planform. The program searches for the minimum 
weight material distribution by searching for the optimum combination of the 
coefficients of the functions describing the material distribution and 
their orientation for fibrous composites. 

F-5A/B CORRELATION?* 

The detailed comparison of skin-thickness distribution, estimated loads 
and structural response used in the WADES design process with those of the 
F-5A/B fighter aircraft was first undertaken. This vehicle was chosen to 
coincide with the ongoing Computer Aided Design Report and Evaluation Study 
(CADRES) currently being conducted by the Advanced Vehicle Concepts Branch 
at NASA/Ames Research Center. It is the intent of this section to validate 
the assumptions concerning the externally applied loads and forces and to 
compare the material distributions and weights obtained in the following 
analysis and design with those of an actual aircraft. In this manner a 
detailed breakdown can be obtained of the factors contributing to discrep- 
ancies in the estimation of pmary structural weight. From this information 
both nonoptimum weight coefficients and areas of analysis or design improve- 
ments are identified. 

In response to a request for structural design information to Northrop 
Corporation, Aircraft Division, a number of pertinent reports were obtained 
from Stanley R. Murnane , Manager, F-5A/B Structural Analysis. These reports 
contained information on the F-5 mass and moment of inertia distributions, 
flutter tests, wing section properties, shear flow and bending stress distri- 
butions, weight and loads data, and tip deflection data. The reports con- 
taining this information are listed in references 1 through 9. 

The F-5A/B correlations are developed generally along the following 
lines: first, comparisons of results from a simplified model with the 
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analyslii of tha actual configuration and ita rasulting rasponaa; and 
arcond, correiationa of the reaulta from a radaaign uaing tha aama geo- 
metric modal and externally applied loada. In the firat part, the F-5 
atructural plan form and aubatructure propertiea were approximated, and 
three critical flight conditiona were aelected. The upper akin-thicknea s 
diatribution waa than aurface fit with the approximate polynomial function 
to be uaed in the WADES progreun to analyze the wing. A aeparate lower 
akin deaign was not conaidered. The program aaaumea equal upper and lower 
akin thicknesses that are computed using an average of the tensile and 
compressive allowable stresses. Actual and approximate wing parameters 
were compared. In the second part, this same configuration was redesigned 
to satisfy the load requirements of the first part. The sjame comparisons 
with the wing parameters of the F-5 were made again. The details and 
comparisons of this procedure follow. 

Summary of F-5 Geometry and Flight Conditions 

Selection of the gt'ometric model and the choice of the critical 
flight conditions for an.xlysis and design of the F-5A/B structure are 
sensitive factors if accurate correlations are to be obtained. Because 
of the current restriction of the WADES program to trapezoidal wing and 
structural planforms, the design tends to be very sensitive to the place- 
ment of the structure itself. Similarly, the selection of the critical 
loading condition directly affects the resulting weight estimate. The 
choice of these loading conditions is often a function of many of the 
parameters in the mission requirements. The particular geometric and 
flight loading conditiona used in this comparison to represent the F-5A/8 
are descrioed here. 

The particular choice of the structural model of the F-5 depends 
upon the positioning of the internal and external configuration of the 
wing. Figure 1 is a cutaway pictorial representation of both the struc- 
tural and non-structural components that make up the F-5 wing. It can be 
seen that the choice of the structural planform is affected by the posi- 
tioning of both the leading-edge and trailing-edge flaps and the volume 
of the wing that ic occupied by the landing gear and aileron operating 
mechanism. 

Figure 2 is the structural idealisation used by Northrop in the 
generation of their internal loads (ref. 2) . Superimposed on that figure 
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also la tha structural planform ussd by th«> *'^AOBS program to modal tha 
P-S. Bacauaa it waa daairad to Includa th# wiftg carry>through atructura 
in tha analysis, tha structural planform was raatrictad to batwaan tha 
15-45-parcant chord linas. Tha nominal axtansion of tha wing tip bayond 
Wing Station (WS) 142.6 was darivad from tha basalina choica of tha F>5 
samispan to bo 151.5 inchaa with tha inclusion of tha Sidawindar (AXM-9B) 
misaila cn tha wing tip. This choic': of the litructural planform waa 
conaidorod to boat modal tha load pat .a of tha major banding loads into 
tha fusalago. 

Tha major diacrapancias of this modal are tha naglact of additional 
structural material aft of the 44-parcant chord line and tha misrapresan- 
tation of the wing-fuselaga junction. The first restriction in modeling 
of the structure represents a 27>percent reduction in equivalent struc- 
tural planform area if only that additional structure outboard of WS 101 
is counted and a 44-perce:;t area reduction if the additional material out- 
board of WS 26 is inclu'jvdi. Tha second restriction at the wing-fuselage 
junction has a twofold effect. It distributes the aerodynamic pressure 
loading over an increase in exposed wing area (27 percent in tha casa of 
the F-5) , and it reacts the resulting shear load at the airplane center- 
line rather than at the wing-body intersection. 

A summary of the F-5A/B wing geometry used in this study is gi^en in 
Tedjle I. !Hie theoretical root chord and semispan were obtained from the 
baseline configuration in reference 10. It is also noted that the F-5 
stores no fi'el in the wing. 

The fligltt conditions used in the correlations with the F-5 were 
chosen to satisfy the critical symmetric maximum wing bending and landing 
loads encountered. Because of internal program restrictions only three 
simultaneous loading conditions derivable from static equilibrium may 
currently be used in a single design sequence. A s’tmmary of the three 
critical flight conditions used is found in Teible II. 

These loading conditions were derived from a combination of infor- 
mation in references 1, 7, and 10. The first two loading conditions reflect 
the identification of the symmetric pull-up and dynamic landing conditions 
in the Group Weight Statement (ref. 10) as being critical. The Wing Stress 
Analysis (tef. 1) identified the first condition as being critical for 
wing stations inboard of WS 114. The second condition was critical in 
sizing of components in the region of the main landing gear trunion. The 
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third condition, a aymmotric pull-up «t lovol, was also idantifiad in 
rafaranca 1 as critical inboard of W8 64 as a rasult of flight tasting, 
nia particular a. angamant of axtamal storas on tha wing for thasa flight 
conditions was primarily obtained from information in tha P-SA/B wing 
Design Loads (raf. 7). 

Rafaranca 1 identified points outboard of W8 85 as being critical 
for several different dynamic store ejection conditions. Tliasa ware not 
included in this analysis because of WADES' inability to reproduce tha 
dynamic loading profile. Tha use of a negative landing gaar weight in 
tha second condition was impleittantad in order to obtain a statically 
equivalent impact load on the landing gear strut. The use of a non-zero 
value for the concentrated loads is Indicative of the positioning of the 
appropriate external store at that wing location. 


Equivalent Cora Properties 


The WADES program does not include a resizing algorithm for the core 
(substructure) properties. The program allows for the input of equivalent 
distributed properties. These equivalent material constants may be obtained 
directly from such materials as aluminum honeycomb or by the calculation of 
a distributed modulus for s spar-rib type of construction. 


Since the F-5A/B used a spar-rib type of internal construction, an 
equivalent density and set of orthotropic moduli were derived to give the 
model the approximate stiffness and weight properties of the actual aircraft. 
To facilitate the development of some averare-distributed properties the 
following assumptions were made; 


e The equivalent material constants to be computed for trapezoidal 

plate in bending with semi-depth, d(^,Ti), are E„, E„, G_, 

X y xz 

G and p core. 
yz ' ^ 

* The relations between actual and distributed core cross sections for 
each of the component moduli are 


Equivalent Distributed 
Cross Section 


Actual y-5 


'Vy'c 


^®^y^ ribs apars®^" ®a 

'®'x> .pars'""®. 
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Equivalent Distributed 
Cross Section 

Actual P-S 

‘“xA'c 

^®^x^ rib webs 

<V\'c 

<®^y) star webs«®**®s 

with 0^.^ "0, and where 
xyc ' 

- sweep angle at x/c ■ 0.3!>, 


• Th« av«rag« F-5 cross-sactional propartlas batwaan WS 64 and WS 89 

ara uaad to computa aquivalant propartiaa in tha y-diraction. 

• Iba avaraga croaa-aactlmal propartlas at x/c « 0.4 ara uaad to 

computa aquivalant propartlas In tha x<>dlractlon. 

• Zn tha calculation of aquivalant bending moduli, both tha flanga and 

wab matarlal In tha spars and ribs are used. 

• Zn the calculation of equivalent transverse shear moduli, only the 

wab material of spars and ribs are used. 

s The aquivalant weight density Is obtained by averaging the estimated 
weight of the spars and ribs over the net volume of the WADES 
structural planform. 

The calculation of the distributed material constants for the core of 
tha P-5 used representative dimensions obtained from reference 1. A 
summary of the estimated equivalent core properties is found In Table ZZZ. 
The details of tha calculations and the values used are found In Appendix A. 

Surface Pit of P-5A/B Thickness Distribution 

To best evaluate the ability of the WADES program to predict both the 
required material distribution and Its corresponding weight In comparison 
with the actual h -5 data, a function was fitted to the actual upper wing 
skln-thlckness distribution. A special-purpose program was written to 
compute the least-squares fit of a ten-term polynomial function used to 
describe the distributed skin thickness. The actual skln-thlckness values 
were obtained for representative locations on the structural planform from 
the summary table of the critical wing loadings In reference 1. The 
details of the calculation of the least-squares fit are found In Appendix B., 
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The resulting functional fit of the actual p-5 akin-thickneea dietri- 
bution is shown in figure 3. Only the basic skin thickness, without the 
additional material due to the presence c*' spar and rib caps, has been 
included in this fit. The contours represent the shape of the resulting 
surface fit. A detailed comparison of the actual and calculated values 
at the input locations is in Appendix B. With this functional fit, a 
complete analysis of the F-5 was performed vithoixt resizing. 

For the previously described structural planform gecmietry and core 
properties, the calculated value of only the structural weight was 507 lbs 
for both wings. This compares to an actual structural weight for the 
F-SA/B of 838 lbs, excluding flap and aileron weights. The calculated 
weight breaks down into '*'^6 lbs for skir material and 221 lbs for core 
(substructure) weight, li.is computed weight results in a 1.65 non-optimum 
weight factor based on the ratio of actual structural weight to computed 
structural weight. If the flaps and other attachments are Included, the 
ratio of the total wing weight to that computed becomes 3.05. 


Since this computed weight was quite low as compared with the actual 
F-5 wing weight, it was further decided to estimate the effect on the 
ccHnputed weights and responses of including the additional material in 
the spar caps 2 md skin-spar attachments. To do this an “effective" skin 
thickness, t^££> was formulated, and a surface fit of this material dis- 
tribution was made. The effective skin thickness was defined to be that 
thickness which would contain the same average cross-sectional area as 
the original skin/spar-cap combination. The equation used to compute 


■eff 


'skin 


At„^^„*w_„ 

caps av 


( 1 ) 


where local p2mel skin thickness, At^^^pg is the nominal 

thickness of the extra material in the spar caps and spar attachments, 
w is the average width of the spar caps, and b is the local width of 
a panel over which the spar cap is to be distributed. A nominal value of 
B 1.0 inches was used in the following analyses. 

The resulting functional fit of the P-5 effective skin thickness 
distribution is shown in figure 4. The contours represent constant values 
of thickness for the resulting ten-term polynomial surface fit. An 
analysis of the F-5 using this functional fit was then made. 

10 




For tha sane geometry and core properties used in the previous fit, 
a structural weight of 552 Jbs was computed for the fit of This 

split into a structural skin weight of 331 lbs and core weight of 221 lbs. 
This is a ratio of actual structural weight to that computed of 1.52 or a 
ratio of 1.69 when compared with the total wing weight (1,041 lbs). 

To obtain a more detailed breakdown of the discrepancies in the weight 
distribution, the spanwise distributions of cross-sectional material area 
and moment of inertia were examined. This was accomplished by chordwise 
integration of material distribution within the structural planform at 
various stations along the span to obtain the structural material area and 
the moment of inertia. To compare directly with the corresponding values 
from the actual P-5, the areas and inertias were resolved into a component 
perpendicular to the x/c ■0.35 reference line. 

Figure 5 is a plot of the computed and actual structural material 
cross-sectional areas versus span. The areas for ''.he actual F-5 were 
obtained from the tidsulated inputs in reference 2. In general, the 
computed values from the functional fits for both t and are below 

the actual values. The major source of error contributing to this iscre- 
pancy is the poor correlation of the actual structural planform and that 
in the WADES model. This is most notice2dt}le at the tip as shown in fig- 
ure 2. An exception is at the root, v4iere the trapezoidal model includes 
more material than necessary. The fourth curve showing the WADES redesigned 
material distribution will be discussed later. Figure 6 is a plot of the 
computed and actual structural moments of inertia versus span. The WADES 
program most noticeedsly overpredicts the inertia at the root where there 
is additional structural planform, euid underpredicts it outboard where 
some planform is excluded. 

Calculation of External Loads 

Computation of the externally applied forces and pressures was broken 
into three phases: the reduction of the given aircraft configuration into 

a set of statically equivalent loads, the estimation of the distribution 
of aerodynaunic forces and pressures, and the summation of these external 
loads into equivalent shear and moment distributions on the wing in a 
form appropriate for comparison. 
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Hie WAOES program accepts as input a breakdown of the aircraft con» 
figuration into the body weight and a set of attached discrete masses, 
and the density of the internal material distribution. Prom this summation 
the program generates the balancing set of aerodynamic pressures required 
for static equilibrium. This static balance of forct^s and pressures is 
graphically shown in figure 7. The aerodynamic pressure distribution, p., 
and the mass distributions corresponding to the weight of the skin, p , 
core, p^, euid fuel, p^, are calculated internally. The presence of a 
landing gear load, input as a negative concentrated mass. (Note 

that stable dynamic eigensolutions cannot be calculated for this flight 
condition.) The center of gravity is assumed to be located at the aero- 
dyneunic center and the balemcing tail load is neglected. Only symmetric 
loading conditions are considered. 

The WADES progreun uses three methods to estimate the external distri- 
bution of aerodynamic forces. These methods are currently available at 
program load time and may not be intermixed. A first-order approximation 
of the loads is obtained by the use of a uniformly distributed constant- 
pressure loading. This loading is independent of Mach number and is the 
simplest to compute. However, it only begins to approximate the loading 
on a thin wing in high supersonic flight. The equivalent pressure loading 
is computed from the ratio of the gross lift required at the maneuver 
loading condition to the wing area: 


„ GLR 

Pa " ~ 

wg 


( 2 ) 


Since the constant-pressure method does not include the flexibility 
of the wing in the loads calculation, a second method for supersonic loads 
generation is available using piston theory to derive the pressure 
distribution. The equation describing this steady-state pressure loading 
is 


P 


a 



Y + 1 
2 


M, 



do _ dw \ 
dx * dx / 


(3) 


where is the angle of attack from the zero-lift line, w is the dis- 

placement of the wing due to flexibility, and y, M^, p^, d, and c are 
respectively the free-stream ratio of specific heats, Mach number, static 
pressure, and the wing depth and ceunber functions. 
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B«cau>« these two aerodynamic loadings are Intended only for super- 
sonic flight, a third method Is needed for subsonic flights. The major 
discrepancy of the first two methods In modeling subsonic flow Involves 
positioning the aerodynamic center at the semi-chord rather than at the 
quarter chord. For the third method, the methods of reference 11 were 
used to derive a modified strip-theory load distribution. The details 
of the corresponding equilibrium conditions associated with this loading 
are developed In Appendix C. 

The static equivalent of a dyneunlc landing condition was derived to 
attempt to model the maximum Izmdlng loads. This load Is approximated from 
a knowledge of the total lift on the wing and the Impact load factor. The 
force on the landing gear then becomes the weight of the aircraft times 
the ultimate load factor less the net lift on the wings at the moment of 
Impact with the ground. Therefore, for both main gear the landing force 
becomes 


F 


ig 


-SGW 



!iift\ 

/ 


+ w 


ig 


(4) 


where ^^Ift fraction of SGW due to aerodynamic lift at Impact. 

The equivalent body weight Is that computed from the statics for the gross 
lift required: 

MCM 

- *‘4“body * “fu.l * «„i„9 + L "cm, * '’l,) 151 

1-1 ^ 

where '"^body only unknown. 

The summation of each of the above external loadings Into equivalent 
shear and moment distributions on the wing was then undertadcen in order to 
present the results in a form for comparison. This summation included not 
only the aerodynamic distributions but also the discrete masses and the 
distributed weight due to skin and core material distributions. The 
resulting spanwise moments and shears have been integrated emd then 
resolved along the 35-percent chord line. This reference line was chosen 
to correspond with the data from Northrop on the P-5. This resolution of 
pressures and forces into shears and moments is typical of beeun modeling 
of wings. Its mesuiingfulness for wings of very low aspect ratio is of 
questionable value for other than a standard of comparison. 
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Die integration of the shears and moments is carried out explicitly 
for the constant-pressure and piston-theory loadings and numerically for 
the modified strip loadings. The integral equations used to compute the 
shears and moments for the constant-pressure and piston-theory loadings 
in the reference axis system along the centerline are 


Yq Ie 

■11 


p dx dy 


p(y - y^)dx dy 


1 TE 


• ■/. /. 


p(x - x(e)ldx dy 


where the pressure, p, is the sum of the distributed pressures due to 
aerodynamic, skin, core, fuel, and concentrated masses, and e is the 
x/c location on the reference chord about which the moments are taken 
(e « 0.35). Similarly, the equations used to numerically integrate the 
shear and moment distributions for the modified strip theory are 




The remaining shears and moments due to material distributions, etc., are 
obtained from explicit integration as in equations (6) . Since the lift 
due to the constant vortex strip is reacted as a discrete load at the 
quarter chord, the resulting torsional moment, My^^, is translated to the 
reference axis as follows; 


”v " ”v ^9 tx(e,y) - X(c/4,y)l 

y y^ Zi 


(8) 


Th« relative senaea of the forcea and moments for the F-5 wirg are shown 
in figure 8. The resolution of the moments into the axis system referred 
to the 35-percent chord reference line is then made by 
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Here the angle d, is the angle of sweep of the reference chord line. 
The calculation and plotting of these distributions has been included 
as an optional output in the WADES progreun. Their calculation is made 
independent of other program functions and does not affect the internal 
force distribution. 


Comparison of Assumed Loads with F-5A/B 

In order to establish the sources of the discrepancies oetween the 
WADES progrsun results and actual F-5 data, a comparison of the calculated 
load distributions with the values of the Northrop wing design loads was 
undertaken. The Northrop design loads used here are summarised in refer- 
ence 7. The original design loads were computed by superposition of rigid 
lift and twist distributions computed for the linear aerodynamic range. 
After completion of the 80-Percent Flight Loads Survey, the measured 
flight data were reduced to provide unit wing shear, moment and torsion 
airload distributions that included the rigid and twist lift distributions. 
These flight data were used for all subsec[uent loads analysis, and the 
original wing-tunnel distributions were discarded. 

Data from only two of the flight conditions input to the WADES program 
to analyze the F-5 are compared here (see Table II) . They correspond to 
the maximum symmetric pull-up at sea level (Northrop #123C-5) and the 
dynamic landind condition (#358B) . A comparison of the spanwise loadings 
computed by the WADES program with constant-pressure loadings and with the 
modified strip analysis is presented. 


Figure 9 shows the comparison of the spanwlse shear and moment dis- 
tributions for the maximum symmetric pull-up case as calculated by the 
WJU>BS program using constant W/S pressure loading, ^o discreparcies are 
noted here. The first major discrepancy is the change in sign of the 
torsional moment, I4y. This is the effect of the assumption of constant 
chordwise pressure distribution. The center of pressure, which is then 
located at the 50-percent chord line, produces a negative torsional moment. 
The second discrepancy is the deviation of the shear and bending moment at 
the root. Where Northrop shows a constant bending moment from the wing- 
fuselage junction inboard, the WADES program shows an increasing value. 

This is due to the failure of the model to account for this interface. A 
secondary effect of this assumption is the inboard t*hi£t of the shear and 
moment curves due to the distribution c£ the pressure loading over the 
entire wing area as opposed to just the exposed wing area. The 16>percent 
increase in pressure corresponding to the difference in theoretical and 
exposed wing area would bring the shear and bending moment much closer 
together. The wrong sign on the torsional moment is unaffected by this 
shift. 

Figure 10 is a comparison of the same symmetric pull-up flight condi- 
tion except that the loads have been computed with the modified strip 
loadings as described in Appendix B. The torsional moment, My, now has the 
proper sense due to the location of the local aerodynamic center at the 
quarter chord. The slight underestimation of the torsional moments is the 
result of improper placement of the chordwise centroids of the attached 
concentrated masses and the distribution of the spanwsie loading over the 
theoretical planform instead of the exposed wing area. A proper choice of 
chordwise centroids would increase the root torsional moment by 100,000 
in- lbs. The recalculation of the spanwise loading, as noted previously, 
would create an outward shift of each of the bending and torsional moment 
curves and of the shear loading. The resulting shear load at the root is 
currently within 5 percent of the actual F-5 data. Only its relative loca- 
tion is in error. 

Figure 11 shows the comparison of the ultimate spcunwise wing loads for 
the dynamic landing flight condition. The results from the calculation of 
loads by the WADES program using modified strip loadings are compared 
with the Northrop loads calculated at the reference time, t “ 162. A 
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landing gaar forca as computad in aquation (4) was usad to astlmate tha 
aquivalant static loads. Tha moments and shaars are generally within 
10 percent of the P>5 data inboard of the landing gaar. However, an 
earlier Northrop reference time, t “ 118, produced the critical landing 
loads. When compared with this case, the current WADES landing loads 
were about 30 percent below the actual dyneunic loads. 

Deflections 

Since the flexibility of the wing almost always affects the distri- 
bution of load over the surface, the spanwise deformation of the actual 
F-5 was compared with that calculated by the WADES program. Figures 12 
and 13 show comparative plots of the spanwise vertical deflection of the 
F-5 wing for a test limit load for Northrop flight conditions #104 and 
#123c-5, respectively. The loads have been reduced to limit load factor, 
and the deflections have been referenced to the aft wing trunion. The 
wing flexibility was calculated using the material distribution obtained 
from the surface fit of the F-5 skin-thickness distribution shown in 
figure 3. 

The deflections at the 35- and 44-percent chord lines are ccxnpared. 
The Northrop spanwise deformations at the 44 -percent chord line were 
obtained from static ground tests. The Northrop deflections at the 35- 
percent chord line are their calculated predictions. In both cases the 
deflections given for the WADES program were estimated from theory using 
a constant -pressure type of loading. It is noted that the Northrop pre- 
dictions always overestimated the actual deformation. The deformations 
calculated by the WADES program using the surface fit of t are within 
2-1/2 inches of the measured Northrop values. The use of the effective 
skin- thickness distribution (not shown in fig.) generally results in a 
10-15-percent reduction in calculated deflections. 


Summary of F-5 Analysis of Surface Pit 
Material Distribution 

Thus far, in order to analyze the F-5A/B, equivalent core properties 
have been estimated, surface fits of the upper skin-thickness and effective 
skin-thickness distributions have been made, and three flight conditions 
have been chosen to model the critical loading conditions. With these 
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models of the actual P-5, the ability of the WADES program to predict the 
structural response and compute the weight of a given configuration was 
evaluated . 

The comparison of the externally applied loads showed that better 
results may be obtained if the modeling of the wing- fuselage junction is 
changed. A switch to a subsonic wing loading from a constant-pressure loading 
produced an improvement in the torsional loading. The deformations computed 
from the surface of the skin thickness showed good correlation with 
Northrop data. The surface fit of the upper skin panels generally was 
within 0.03 inch of the actual skin thickness at any point on the wing. 

The integrated cross-sectional area distribution perpendicular to the 
35-percent chord line for the effective skin-thickness function showed 
the beet results inboard of the landing gear. However, poor correlation 
was obtained outboard of that spanwise station due to improper modeling 
of the structural pl 2 uiform. The total weight computed from a surface fit 
of the actual F-5 skin was significantly in error becf.iuse of the presence 
of “non-optimum** weight and because of the reduced structural planform. 

Redesign of F-5 Thickness Distribution 

In the first phase of the F-5A/B study, each facet of the structural 
and aei.odynamic analyses used by the WADES program to model the structure 
and loads was compared. In order to evaluate the design capability of 
the progreun, a redesign of the original thickness function was undertaken. 
This entailed designing the wing for strength using the thickness function 
and loads of the previous analysis as the starting point. The previous 
geometric representation, flight conditions, emd equivalent core properties 
were used. During this design the coefficients of the function describing 
the thickness distribution were optimized to obtain the minimum-weight 
structure to satisfy the strength and m ..nimum-gage constraints in the wing. 

Though the F-5 wing presum2d>ly satisfies all the design requirements, 
because of modeling differences such as using a distributed core and 
thickness function and a different structural planfoxrm, the analysis model 
with the initial surface fit did not satisfy the set of WADES constraints. 
This is apparent in figures 9 auid 10, where the bending moment at the root 
significantly exceeds that of the actual F-5 due to improper modeling of 
the wing-fuselage junction. The resulting redesigned thickness distribution 
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ov«r«stiroat«s tha matarlal raqulramants at tha root. The plot of tha WM)BS 
radaalgnad croaa-aactlonal araa of atructural matarial varaus apan in 
flgura 5 axcaada tha actual araa at tha root but at ill undaraatimatoa it 
at tha tip. 

Tha raaultlng radaalgnad P-5 wing waighad 670 Iba whan daaignad for 
the conatant-praaaura loadlnga. Of that walght, 449 Iba was in tha upper 
and lowar skina. Tha cora weight waa tha sane as before. That total 
computed weight yields a non-optimum weight factor of 1.25, based on tha 
ratio of actual structural weight to computed atructural weight. A 
similar design using the loads computed from subsonic strip theory (in 
fig. 10) weighed 602 lbs, reflecting the change in spanwiso load distri- 
bution. Though these weights appear to be closer to the actual P-5 
weights, their spanwisa distribution is actually worse than the surface-fit 
values. The major part of this discrepancy comes from the difference in 
bending mcMnents at the root as evidenced in figures 9 and 10. The displace- 
ments of the redesigned thickness distribution are edsout 8 percent less 
than the effective skin-thickness displacements. 

GENERAL WEIGHT CORRELATION 

Por the WADES program to provide accurate estimates of wing weights 
over a broad class of aircraft, a general non-optimum weight factor based 
on the ratio of actual wing weights to computed values must be statisti- 
cally estzdslished. Previously, such a ratio was computed for a specific 
aircraft, the P-5A/B. In this section such a factor will be computed for 
a class of aircraft and a regression analysis performed to attempt to 
establish its value. 

The class of aircraft used in this wing-weight correlation consists 
of U. S. Air Force and U. S. Navy fighters. This grouping was chosen 
because of the adaptability of medium-to-low-aspect-ratio aircraft to the 
plate theory structural model used in the analysis. Because of the high 
performance requirements of these aircraft, their designs display a high 
dependence on the strength requirements and loads imposed on the wing 
structure. Because mamy were boarderline on incurring weight penalties 
for required aeroelastic stiffening, the computed "non-optimum" factor 
may also reflect such additional material. 
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Th« approach takan to davalop this walght corralation factor waa to 
first parform a daslgn of tha givan configuration baaad on tha bast astl- 
mats of tha structural modal, and than sacondly to parform a llnaar 
ragrassion analysis with tha computad optimum walght to find tha bast 
factor or factors to corralata with tha axlstlng aircraft wing %^lghts. 
Savaral combinations of wing component weights ware triad to asse<‘'°i tha 
dapanda.icy of tha total wing walght on them. 

Preliminary examination of tha estimated component and total wing 
weights Indicated that tha bast correlations ware obtained by comparing 
only tha weights of tha structural planform. Because tha function 
describing tha thickness distribution Is continuous over tha entire plan- 
form, tha computad weights are misleading In that they do not account for 
the discontinuity In material between primary and secondary structure such 
as flaps and ailerons. As a result, tha Integrated material volume Is 
grossly overestimated In these regions. Subsequent correlations ware made 
using only the weights of the material contained In the structural planform. 

Aircraft Used In Weight Regression 

The aircraft < sed In the wlng-welght regression analysis tc determine 
the non-optimum weight factors were U. S. Air Force and U. S. Navy 
fighters. Tables IV emd V contain a summary of the wing parameters for 
the vehicles considered in this analysis. Table IV containc a list of the thick- 
ness-to-chord ratio at the wing root and tip, the root chord, semispan, 
leading-edge and tralllng-edge angles, an approximate chord fraction of 
the leading-edge and traili ig-edge structure, and an estimate of the frac- 
tion of the structural planform available to contain fuel for each of the 
aircraft studied. Because the structural planform does not always align 
itself with the constant chord lines assumed by the program, chord fraction 
of leading-edge and trailing >edge structure was selected to approximate an 
eqxiivalent structural plemform ?~rea. If the volume fraction of the avail- 
able fuel was not known, a default value of 0.5 was used, if the weight 
of the fuel in the wings Is known, this value may be later computed. 

Table V Is a summary of the critical loading conditions input to the 
WADES program to design the various aircraft wings during the correlation. 
The component weights and load factors were obtained from the vehicle 
group weight statements in reference 10. Th<i F-5 loading conditions are 
a composite of the group weight statement data emd information in 
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r«f«rMic« 7. Wh«r« th« design altltuda and Mach nuinbar wara not availabla, 
dafault valuaa of 25,000 faat and 0.9 wara uaad. Tha uaa of nagativa con- 
cantratad loada indlcataa application of tha landing whaal load aa computad 
in aquation (4) . 


Ragrasaion Analysia 

Bacauaa of a lack of diract corraspondanca batwaan computad and actual 
wing waighta during praliminary axaminat.'on of dasign waights, four sapa- 
rata ragraaaion analyaaa wara mada to datermina tha bast non-optimum 
waight factors basad on various componantn of tha structural waight. Zn 
addition, tha waights wars computad for two typas of appliad loadings: a 

constant -prassura load, and tha subsonic modifiad strip loading dascribad 
in Appandlx C. *n\a corralatlon factors ganaratad hara wara obtainad on 
tha basis of tha total wing waight, including tha additional control -surface 
weight. Another valid non-optimum waight constant might bo computad basad 
only on tha ratio of actual to computad structural waight. This was not 
undertaken hara, since tha astimation ot tha total wing weight was of 
primary Interest. 

The variables used in this regression emalysis ware the waight of the 
skin structure, WTSS, the weight of the core structure, WTCS, and the plari- 
form area of the wing, The weight of the core structure was computed 

from the product of the average core density and the volume of the core. 

S^,g waj included to ch^ck for the dominance of control surfaces and sub- 
structures. The four equations used to fit tha wing weight data were; 



- B(1)WTSS + B(2)WTCS 

(10) 


- B2 WTSS WTCS 

(11) 


- B3 WTSS 

(12) 

w 

4 

- BS(1)WTSS •«- BS(2)S /268 

wg 

(13) 


All weights are in pounds and is in square inches, in all cases the 

weight of the skin structure was used as one of the independent variables. 
Percent errors based on the actual wing weight and the weight computed in 
equations (10) to (13) were computed as follows: 


P 


- 1), i - 1,2, 3,4 


(14) 
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r«gr«asion analysis for sach of tha abova aquations was parfortnad 
usinq a laast-squaras functional fit. 

nia four aquations usad in tha ragrassion analysis wars salactad in 
ordar to avaluata tha sansitivity of tha computad waight to tha various 
componants. Tha waight of tha skin structura was usad in aach of tha 
aquations, sinca it contains tha only componant dasignad by tha WAOBS 
program and carrias tha primary load in tha wing. Tha first aquation was 
salactad to astablish tha ralativa significanca of tha waight of tha cora 
(substructura) . Bacausa tha cora voluma is diractly proportional to tha 
voluma containad within tha structural planform, tha fraa coafficiant on 
WTCS than bacomas an astimata of avaraga dansity of tha subatructura. 
lha ragrassion analysis using tha sacond aquation was mada on tha basis 
that tha valua of tha dansity of tha core was computed from tha astimation 
of tha aquivalant distributed properties of the P-5 as derived in Appendix A. 
The third equation was used to evaluate whether the wing waight was diractly 
proportional only to tha weight of the skin. Since most minimum-gaga 
effects in non-primary structura are proportional to the planforro area of 
the wing, the regression analysis using the fourth equation was performed. 

The parzunetars used to describe the geometry of the wings were not 
included in the regression analysis. It was assvimed that their effect was 
included implicitly in the design of the wing itself. 

Regression Results for Isotropic Wing Design 

The regression analysis using the four weight equations was carried 
out for the two types of static loading discussed previously. In each case 
the minimization of the weight of the skin structure was taken as the objec- 
tive of the design. The results of the regression analysis for the weights 
computed using constant-pressure loads are shown in Table VI. The results 
generally show a dominance of the core weights and the term proportional 
to the planform area. The third equation demonstrated the weakest 
correlation. Table VII contains the results of the regression emalysis for 
the wing weights computed using the modified strip loads. The results show 
a strong dominance of the weight of the skin structure in the first three 
equations. The wing area still exhibits a strong correlation in the fourth 
expression. 
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Upon consldoratlon of th« weights computod from th« ourfaco fit of 
tha thlcknaaa dlatributlon of tha F-5 wing and th* .^clanciaa In modallng 
tha atructural planform, it ia racommandad that tha firat rograasion aqua- 
tion ba uaad with tha coafficianta B(l)"2.6 and B(2)"0.7d. Thaaa valuaa 
corraapond to coafficianta computad for tha waighta aatimatad uaing tha 
modifiad atrip loada. lhay alao compara vary cloaaly to tha valuaa that 
would ba coiqputad from tha aurfaca fit of tha affactiva akin-thicknaaa 
diatribution. lha currant program daficianniaa in modal ing tha atructural 
planform auggaat that tha aatimataa of tha waight coroponant proportional 
to cora voluma may ba arratic; and tharafora, tha lowar corralation 
oafficiant for B2 ahould ba uaad. 

Ragraaaion Raaulta for Compoaita 
Matarial wing Daaign 

An attampt waa made to obtain datailad information on apacific wing 
cMRF'^nenta built from cong>oaita material. Tha information racaivad did 
not contain aufficiant information to check waight aatimation directly. 

In lieu of apacific data, an alternate procedure baaed on obtaining a 
fixed percentage reduction in weight over axiating aircraft waa undertaken. 
Weight correlation factora baaed on the apacified percent reduction ware 
then obtained from the previoualy deacribed regreaaion analyaia. 

To iti.plement thia procedure the group of U. S. Air Force and Navy 
fighter aircraft wereredeaigned uaing compoaite materiala in the wing 
cover panela. The atructural akir and core weighta were recomputed. No 
modification of the eatimated core density waa made. The current wing 
weight waa multiplied by a constant fraction, and the regression analyses 
using equations (10) to (13) were carried out. 

The correlation coefficients were determined for three wing weight 
percentages: 100, 70, and 60 percent of the original wing weight. The 

wings were designed using constant-pressure loadings. The results of the 
regression analysis for an eatimated 7— percent wing weight are shown in 
Table Vlil. The regression results for 100 and 60 percent were generally 
within a constant of these values. The coefficients shown here do show a 
strong correlation with the computed akin weights. If the correlation 
coefficients developed for the design of wings with isotropic wing skins 
are used, the weights are generally between 60 and 70 percent of the actual 
vehicle weights. In view of the current state of technology, it is 
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r«c«nmondf)d that the slightly higher coefficients derived here should be 

tised . 


INTEGRATION OP WAGES PROGRAM INTO ACSYNT 

The primary programming task accomplished in this study was the adapta- 
tion and incorporation of the prior work of Dr. Mullen, involving structural 
optimization techniques for automatic resizing of low-aspect -ratio wings, 
into Ames Research Center's ACSYNT program for vehicle synthesis. This 
task was completed and a teat case involving the convergence of an aircraft 
to satisfy the mission requirements of the P-5A with ACSYNT was run. Some 
of the modifications required to interface the WAGES program with ACCYNT 
are summarized here. 


Operation in ACSYNT Overlay Structure 

The primary programming task to Integrate the WADES program into 
ACSYNT v;as its conversion to an OVERLAY structure. This was necessary in 
order for the WAGES program to reside simultaneously in core with ACSYNT 
within the CDC 7600 core limitations. The WADES program was sufficiently 
modularized so that no major adjustment in the progrzun flow was required. 

The problem then became one of maintaining as many of the program features 
as possible without sacrificing progreun generality. 

In arriving at the current overlay structure two subdivisions of the 
program were considered. In each case the WAGES program was required to exist 
as OVERLAY 5,0 within the ACSYNT overlay structure and maintain its own sub- 
overlays. In the first attempt, the program was set up with an executive 
main overlay to branch to the appropriate function depending on the 
request from ACSYNT, and three sub-overlays which provided three functions: 
input, analysis or design, and detailed output. This information was 
provided according to the requ.est for information from the ACSYNT param- 
eter ICALC. This breakdown provided the most direct program flow with 
the least exch^mge of overlays in and out of the machine during execution. 
This version was made operational initially for only the strength design 
of isotropic wings. Because of the heavy demand for core space at that 
time for the analysis and design overlay, it was determined that it would 
not be possible to have the composite strength design code or its stiffness 
and flutter design code reside in core without significant reduction in 
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th« cora storage required by ACSYNT overhead. The initial load with this 
version required aOsout 168K octal core locations to load. This limit 
proved to be unacceptable in view of machine access requirements imposed 
by the computer operating system. 

The second and current overlay structure again used OVERLAY 5,C mainly 
as an executive function with the sub-overlays providing the three func- 
tions: input/output, wing analysis, and design for strength only. In 

this manner the analysis code and the design code, which used CONMIN 
(program for CONstrained function MINimization) , could be separated. The 
main routine that organizes the sequence of analyses and designs also 
resides in the same overlay as the WADES executive routine. *nie loss of 
generality resulting from this choice of overlays is the inability to 
perform a combined strength and flutter design simultaneously with the 
ACSYNT program. A flutter analysis is still possible in this mode; however, 
the rec[Uirement that the analysis and design code both reside in core pre- 
cluded this method of operation. The possibility of design for flutter is 
still available in a stand-alone mode. The decision to implement this 
overlay breakdown was made on the basis that the only mode of operation in 
which the WADES program would be used with ACSYNT in the near future would 
be in the design for strength only. A return to the first overlay structure 
outlined above, to permit combined strength and flutter design, would be 
possible given a 25K octal reduction in OVERLAY 0,0 core requirements. 

The flow of calculations through the wing design executive routine, 
OVERLAY 5,0, is controlled by two parameters, ICALC and ICCftITR. The first 
is the ACSYNT control parameter, and the second is a user-specified 
control input. ICONTR determines the branching to either an analysis-only 
mode or to the wing design for either an isotropic or a multi-layered 
composite wing. The flow chart in figure 14 outlines the basic subroutine 
and OVERLAY flow of the WADES executive routine STRUCW (OVERLAY 5,0) with 
the branching determined by the two control parameters, ICALC and ICONTR. 
Similarly, in figures 15, 16, and 17 are the basic flow of OVERLAYS 5,1, 

5,2, and 5,3, respectively, and their corresponding subroutine usage. 
Included in figure 15 is the branching according to the value of ICALC, 
and similarly in figure 17 is the branching as determined by ICCXITR. In 
figure 16, the basic subroutine usage in the analysis of the wing stiff- 
ness, loads, stresses, etc., is shown. The call to the various analysis 
routines are determined by the control parameter IANAL(I,IFLT) as required 
by the IPLT'th flight condition. 
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The design of both the isotropic and multilayered composite wings 
employs a sequence of minimization problems to arrive at the design weight 
of the wing. The relative error in the skin weight below v^ich this con- 
vergence is forced is at the discretira of the user. lAie basic subroutine 
and overlay flow of the WADES program as it exits with ACSYMT is shown in 
figure 18 for the design of isotropic wings and in figure 19 for the design 
of multi-layered composite wings. The number of iterations necessary to 
obtain the wing weight may be specified either as the nunOder of iterations 
required to converge the weight to within the desired error or as a 
maximum number of iterations (NRAT) . 

Geometry Interface with ACSYMT 

In order to interface the geometric descriptors of the WADES program 
with those of ACSYMT, the equations describing the relationships among the 
appropriate variables were derived. The W/>DES program is currently 
limited to tr?ipezoidal planforms and uses the root chord (R) , semispan 
(SPAN) , and leading-edge (THETl) and trailing-edge (1HET2) emgles to 
describe the geometry. ACSYMT, on the other hand, uses a nondimensionalized 
description with an arbitrary reference line, which is usually the quarter 
chord. The remaining variables, t/c at the root and at the tip, are 
identical. 

The basic external planfcrm description of the wing used in the WADES 
program is shown in figure 20. The equivalent gecxnetric values were derived 
in terms of the appropriate ACSYMT descriptors. The WADES geometric values 
are summarized in terms of their ACSYMT equivalents in figure 20. 

Loads Interface 

In order to interface the loading conditions used by the WADES program 
with the changing weight and flight information generated by ACSYMT, some 
method had to be devised which could update the loads during execution. 
Further, this method had to reflect the nature of the critical structural 
design conditions and not necessarily just the mission flight profile. 
Because of program limitations, it also had to be limited to a maximum of 
three such critical conditions. Three approaches are outlined here that 
may be taken to generate the n.'cessary structural design conditions to be 
met. 
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The first approach is to read in the critical structural flight con- 
ditions directly. This method was used in generating the correlation 
information presented later in this report. In this case the atmospheric 
and weight information is included in the input data. In the particular 
instance of the vehicle correlations, these design conditions were obtained 
directly from the Group Weight Statement furnished by the airframe manu- 
facturers. These included the structural gross weight, ultimate load tac- 
tor, and fuel contained in the wings at the design conditions. The ccndi- 
tions generally contained the design flight and landing weighs, maximum 
gross weight with zero fuel in the wing, a catapulting condition, where 
appropriate, and the minimum flying weight. Usually, the particular arrange- 
ment of external stores and maneuver conditions for the configuration were 
not included. In most cases, an estimate determined by working back from 
the gross weight and fuel condition to a configuration had to be used. 

This approach of reading in the design flight conditions is valid 
only if a known configuration is being analyzed. Even then, it remains 
true only if the remaining body and fuel conditions are constant through- 
out; the design. Since the convergence portion of the ACSYNT progrzun 
operates in a mode where most of the individual components are continually 
being updated, a direct input of the structural docign flight conditions 
would lead to erroneous results. 

The second approach to providing the critical structural flight 
conditions is to incorporate a special-purix^se subroutine to compute the 
appropriate loading conditions. In this manner the structural loads czm 
be made to reflect the particular configuration and mission requirements. 
Such loading conditions can then be made to include such items as partial 
fuel conditions in the wing and the positioning of empty fuel tanks on 
the wing during lemding. Ihis method is used later to update the loads 
for the F-5A/B case study, checking the convergence of the WADES program 
with ACSYNT. 

In the F-5A/B study a specific routine , FLTLDS, was written to 
specify the approximate flight conditions outlined in the Group Weight 
Statement (ref. 10) . Instead of specifying the loading condition directly 
from the design structural flight conditions, the critical loads were 
written in terms of the general gross takeoff weight and various fuel 
weights. Only the relative proportions of fuel at maneuver and Izmding 




27 


were maintained similar to those in the weight statement. Thus, the 
critical stress gross weights and loads were updated during the ACSYNT 
convergence cycle. The particular relationships used to compute the 
appropriate WADES loads inputs are summarized in Tiable DC. Though the 
maximum symmetric pull-up and the landing load condition are two of the 
major design considerations, routine PLTLDS is P-SA/B aircraft-dependent 
because of the inclusion of such factors as the values of the fuel 
fractions considered and the positioning of fuel tanks on the wing tips. 

The third approach to the specification of the structural flight 
conditions would be the derivation of a general routine to find the 
critical conditions. In general this would entail a survey of the 
critical gust, maneuver, landing, and eventually flutter conditions. It 
would also require a check of the possible external store configurations. 

An interim approach would be to utilize the outline of a routine such 
as the F-5A/B -derived PLTLDS with the ability to input all fractional 
relationships. If wing-mounted engines are to be considered., their 
positioning on the wing should be included. At the moment no such geo- 
metry descriptor exists within ACSYNT to locate their chordwise and span- 
wise locations on the wing. Such interfaces would have to be either 
generated or input to obtain their impact on wing design. 

Optional Material Properties 

The WADES progr 2 un uses two modes to input material properties for use 
in analysis or design. In the first mode, the material properties are 
read as part of the normal input data stream. In this case the elastic 
constants and density are input through the namelist MATERL and the failure 
stresses are input through the nzunelist CNSTR. in the second mode, the 
WADES program generates the required properties and failure criteria 
internally. In this case the appropriate material constants are defined 
for three materials: (1) aluminum, (2) titanium, and (3) graphite/epoxy. 

The last defines only the appropriate lamina properties and is used in 
multilayered composite analysis and design. 

The second mode of material property input is availed^le optionally 
by the input of a non-zero value of the prograun parameter, ITYPES, in 
namelist OPms. In this mode of operation the appropriate material 
properties are defined prior to their input in the normal data stream. 

Thus, any of the isotropic material constants defined in this manner may 
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btt overwritten by Ita appropriate redefinition in the namelists MATERL 
or CNSTR. The definition of these material properties is performed in 
subroutine MATRLS. A summary of the values available by option is shown 
in Table X. 


Special Version of CONMIN 

The WADES program uses the mathematical -programming technique called 
the Method of Feasible Directions (MFD) to perform the optimal resizing of 
the wing structure. This technique was originally outlined by Zoutendijk 
(ref 12) and demonstrated for use in structural optimization in refer- 
ence 13. This MFD algorithm has been programmed in a general form in 
the Fortran progrzun for CONstrained function MINimization (COMMIN ) , 
reference 14. 

Since this program also acts as the resizing algorithm for aircraft 
optimal design in the ACSYNT control program, it was necessary to include 
a second version. This version provides a sub-level optimization function 
and has to reside simultaneously in core. To avoid a Fortran naming 
conflict the name was changed to CONMM with subroutines CNMNJl, . . . . ,CNMNJ9. 
This version was redimensioned to handle thirty design variables and up 
to forty active constraints. 

Convergence of F-5A/B with ACSYNT 

To check the operation of the WADES program with ACSYNT the two 
programs were connected and several test cases were run to test their 
convergence properties in an iterative mode. The connection to the October 
1974 version of ACSYNT using a temporary buffer routine and the F-5A/B 
version of FLTLDS was made. In figure 21 are the results of a WADES/ACSYNT 
convergence cycle. In the figure the convergence characteristics of three 
typical variables, WWING, WFUEL, and WGTO, are plotted versus the iteration 
number. The wing weight, WWING, in this figure was computed by the WADES 
module. Tlte remaining two variables, the fuel weight and the gross tadceoff 
weight, were computed in the current trajectory and weights estimation parts 
of ADSYNT. This wing weight was computed using the preliminary correlation 
relationship 

WWING - 2.4 WTSS + WTCS (15) 

where WTSS is the weight of the skin structure and WTCS is the weight 
of th'i internal substructure as c^puted by WADES. These correlation 
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coefflclants do represent a good value for the F-SA/B using constant W/S 
wing loadings because of the additional detail comparison, but they are 
not the set of correlation coefficients as determined earlier in this 
final report for a broader class of aircraft. 

A detailed examination of the convergence cycles shows that the 
trends exhibited by the wing weight do follow closely the weight of fuel 
and gross takeoff weight. The first five iterations in figure 21 are 
spent bounding the limits of WTSUM-WGTO. In the next three iterations 
the design has essentially converged. A total of twelve vehicle analyses 
were used to converge the aircraft to the necessary tolerance specified 
by ACSYNT. An average of four analysis and design cycles were used by 
WADES to converge a wing design during each of the ACSYMT iterations. 

This average convergence rate should drop significantly when redesign, 
which uses the previous design for the starting point, is iropleinented . 

The horizontal line on each curve in figure 21 represents the initial 
estimate of each particular parameter, which is equal to the actual value 
obtained from the F-5A/B Weight Statement. In each case the final weight 
does converge to a value lower than the actual. It is noted that the wing 
weight was slightly overestimated on the first iteration when this set of 
correlation coefficients was used. 

WADES PROGRAM DOCUMENTATION 

The purpose of this section is to describe the basic usage of the 
wing aeroelastic design program, WADES, with the ACSYNT progrcun for 
vehicle synthesis and in a stand-alone mode. Included are the basic 
description of parameters, a sample input, a sample of the printout from 
routine WOUT, a short description of the purpose of each routine, and a 
description of the use of alternate analysis and design routines. The 
program computes the stiffness e.nd mass properties for a wing using an 
"equivalent -pi ate" Rayleigh-Ritz model. The structural response is cal- 
culated for the application of both steady and unsteady aerodynzunic loading, 
and the material distribution of the wing skin cam be resized to satisfy 
both strength and aeroelastic requirements. The theoretical background 
for the development of the analytical model was originally determined under 
sponsorship of the National Research Council, and the report summarizing 
this effort is in preparation. 
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Description of Parameters 


The geometric description, material properties definitions, flight 
conditions, and design constraints are detailed in this section. The 
descriptions in this section pertain to those definitions affecting the 
determination of input parameters. All information required to derive the 
inputs for both the stand-alone progrzun version and the integrated version 
used as a module for the ADSYNT program is explained. Only minor modifi- 
cations in several routines have been made to adapt the WADES progrzun to 
operate in the ADSYNT overlay structure. The current version of the stand- 
alone program version will execute with all options in 142K octal words of 
core. 

Since both programs use essentially the same routines, only minor 
omissions in the input data must be made to execute the WADES program with 
ADSYNT. Both programs use the same input subroutine. Only the variables 
describing the geometric shape of the wing and those weights which vary 
with changing gross weight need to be omitted. In subsequent analyses 
those values v’ill be overwritten by values supplied by the ADSYNT main 
program. A basic user's guide to the WADES program inputs is given in 
figure 22. 

Wing geometry .- The geometric planform analyzed in the WADES program 
is trapezoidal and consists of superimposed aerodynamic and structural 
regions. The structural planform is always contained within the aerodynamic 
planform. Ail material contributing to bending strength is contained within 
the structural planform. Nonstructural material within the aerodynamic 
planform is considered to contribute only to wing weight and mass properties. 
Core properties are considered to be distributed over both structural and 
aerodynamic planforms. Fuel is considered to be distributed only within 
the structural planform. 



where 

R « root chord (length) 

SPAN • aemispan (length) 

THETl - leading-edge sweep (degrees) 

1HET2 ■ trail ing-edge sweep (degrees) 

XLE - location of leading edge of structural planform 
(fraction of chord) 

XTE - location of trailing edge of structural planform 
(fraction of chord) 

Thickness and depth functions .- Thickness, and depth and camber distri- 
butions are represented as continuous functions. The particular depth or 
Ccunber function used in the program is a polynomial with zero depth enforced 
at the leading aid trailing edges. This results in a symmetric airfoil 
section with sharp leading and trailing edges. The resultant function is 
formed as the product of the planform polynomial, WP(|,ti) , which enforces 
zero depth at the edges, and a user-supplied polynomial, PD(|,ti). An 
approximation to a biconvex wing section is available as a default within 
the program. The shape of the ceunbered surface is also specified in the 
same form as the depth function through the variable, PC(^,t)). The thick- 
ness is similarly the product of the planform polynomial and a polynomial 
distribution function plus a minimum thickness constant. For multilayered 
composite design each lamina may be described by a separate function. 

These polynomial functions may then be written: 
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Planform polynomial - WP: 

WP(^ . n) ■ -n*TAN (THETl) *SPAN/RN) * (-^+t)*TAN (THET2) *SPAN/RN+R/RN) 


Depth and thickness functions - D,T: 

D(^,ti) - WP(?,n)*(PD^ + + ... + pd^pd?^"^ 

+ PDgTi + PD^^g^T] 

s 

+ pd^t)^“M 

i-i,IPD; j-x,jPD 


T(^,n) - WP(?,ti)*(PT^ + PTi^i^ + ... + 

+ PT n + FT . ^ f n 
. 2 ‘ 1+2^ ' 


+ PT.t)^"M + T 

1 I 


min 


i-x ,IPT; j-x ,JPT 


Alternate depth 
tip 

XTCR 

XTCT 


representation (approximate) : 

thickness to chord ratio at root 

thickness to chord ratio at tip 

location of maximxiro t/c along root 
(fraction of chord) 

location of maximum t/c along tip 
(fraction of chord) 


I 


This approximates a linear t/c and linear XTC distribution to the 
surface of the wing as a depth function. The number of degrees of freedom 
it this approximation is determined by IFD. If XTCR"XTCT«0, a biconvex 
section is fit with JFD-1 and NFD>IFD. If XTCR>0 , a linear fit to the 
location of maximum t/c is included with JFD»2 and NFD"2*1FD-1<11 . 

Material property definition .- All material properties are considered 
to be distributed and continuous throughout the wing section. All properties 
are given with respect to the global reference system. Isotropic material 
properties are considered constant throughout the planform. Composite 
properties similarly maintain the- properties of the discrete ply over the 
entire planform. Core properties are also considered as distributed and 
continuous throughout the section. The core may be modeled as either 
conventional or sandwich construction and appropriate properties averaged 
through the section. 



Material properties 

Skin: EXS ,EYS ,GXYS ,GXZS ,GYZS ,RHOS 
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Loads . - The effects of inertial, discrete, and distributed pressure 
loadings are considered. The discrete loads are incorporated as concen 
trated masses. Their effect is included as a discrete inertial force 
loading in the static analysis. They are included as discrete inertial 
masses in the generation of mass properties. Discrete loads are located 
at fractions of chord and span. Allowance is also made for the inclusion 
of certain weights (such as external fuel or armament) as a function of 
flight condition. Similarly., the effects of the inertial loading of the 
skin, core, and fuel weights are included in both the static analyses and 
mass properties. 

Inertial loading of discrete loads: 


PCM(^,n) - -ANZ Z WCM(XCM,YCM )6 (XCM) ,YCM) 

I-1,N04 


3 


Similarly as a function of flight condition for 


J-1,N07C 

Z-1,NCMPLT 


martial loading of skin, cora, and fual: 

P, (C.n) • -AWZ P,T(^,T|) 

ovar aarodynamic pi an form 

Pc<5.n) - -2 ANZ PcD(^.n) 

ovar aarodynamic planform 

P^(^.n) ■ -2 ANZ pjPVA D(^,t)) 

ovar structural planform 

where 


ANZ - ultimate load factor 
p . • densities of skin, core, and fuel (#/L**3) 

PVA ■ fraction of volume available containing fueKPVOL 
WCM ■ discrete inertial load 
T(^,r|) ■ thickness distribution of skin cover plates 
D(^,t]) " depth distribution of core 

Distributed pressure loads are incorporated in two forms: second- 

order piston theory and constant-pressure loading. The piston theory 
pressure distribution is computed iteratively as a function of the angle 
of attack required for gross lift and the displacement shape. A constant- 
pressure distribution is availedsle as an alternate loads subroutine. This 
routine computes a constant-pressure loading from W/S. 


Pa^»l) - GLR/SWG - W/S 


S^cond^ordar piston thsory prsssu^*s distribution: 



Ths tsrm including ths rats of changs of displacsmsnt with tlms prov/.dss 
ths basis for ths asrodynamic forcing or damping ia ths calculation of ths 
fluttsr Mach nui^sr and frsqusncy conditions. 

A modifisd strip thsory for subsonic, static loads is also availabls. 
Thsss loadings par unit span urs intsgratsd and applisd as discrsts forcss 
along ths quartsr chord. This form of ths asrodynamic loads is bassd on 
ths msthod of Orsy and Schank (NACA ISl-SOSO , ref. 11) and is restricted to 
small angles of attack at subsonic speeds. 

Ths general msthod for computing ths static balance of external 
forces is by the satisfaction of the requirement that the gross lift 
available, GLA, is equal to the gross lift required, OLR. Two methods 
of computing the gross lift required are included which are independent 
of the type of pressure loading used. The choice of method of calculation 
is specified by the parameter, IAMAL(2,IPLT) , for the given flight condi- 
tion. For XANAL(2, IFLT)"1, OLR is equal to the summation of body, skin, 
fuel and concentrated maas weights times the memeuver load factor. 

OLR • AMZ(WBODY-Hif'lVING+WFR-fNWINGS(SUM WCM(I) 

+ SUM AMFLT(I) ) 

For IAMAL(2,1FLT)*2, OLR is computed from the input of the specified 
wing loading required. 

OLR - (W/S) *SWG 

The WADES program is organized to consider NOFC flight 
environments and loading conditions during both analysis and design. 

The speed of sound, Mach number, static pressure, and ratio of specific 
heats are the environmental factors which vary with flight condition. 

The discrete loads, load factor, body weight including the payload and 
fuel not in the wings, and fuel weight in the wing are the loading 


37 


conditions which vary with flight conditions. The encl-ssed program version 
is dimensioned to handle three such flight conditions of various combina- 
tions of the above environmental and loading conditions. 

Behavioral constraints . - nie behavioral constraints considered are: 
stress, displacement, fundamental frequency, and flutter Mach number, 
frequency, and dynamic pressure. The behavioral constraints are checked 
at each required flight condition. For an isotropic material the 
strength constraints are the maximum allowable von Mlses* stress resultant, 
SMAXT. The stress is evaluated at a grid of MXSIG by NYSIG points over 
the wing planform. This grid is restricted to points within the structural 
planform of the wing. The displacements of both the leading and trailing 
edges at the wing tip are constrained to be less than the allowable, WMAX. 
The fuudeunental frequency is constrained to be above its minimum, EIGMIN. 
The flutter frequency, Mach number, zmd dynamic pressure must be greater 
than their corresponding minima, FFMIN, FMMIN, and QFMIN. 

For multilayered composite materials, the strength constraints are 
in the form of a modified distortional energy criterion for the failure 
of the individual plys according to their direction of orientation. This 
constraint is evaluated at each of the points in the structural grid. This 
constraint teUces the form 


1 > 



-KL12 


SLll 



wher_ the SLIJ ar* the uni -directional failure stresses for tension or 
compression in the appropriate direction (ref. 15) . 

Geometric constraints .- The geometric constraints considered are for 
minimum gage material thickness of the wing cover sheets. The thickness 
function is written as a function plus a constant. The constant has been 
prescribed as the minimum gage and the remaining functional is constrained 
to be positive at all points. The functional is evaluated at a grid of 
NXGC chordwise by NYGC spanwise points over the planform of the wing. 
Th'S number of points in each direction must be greater than the order of 
the polynomial in the appropriate direction to insure a non-trivial 
solution for intermediate points. 
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For multilayered composite materials an additional constraint of 
the material stiffness in the various principal directions of the composite 
has also been imposed because of the highly orthotropic nature of the 
fibrous composites. This constraint takes the form 


N P L Y S 

^ml J J 

where t^^ is the thickness of the ply and ^£iij is the consti- 

tutive relation of the lamina in the i^^ component direction of the 
composite for the ply* 

For che design of isotropic material wings for strength, only the 
von Mises* stress and minimum gage constraints are used. Similarly, 
for the design of leuninated composite wings, only the modified distortional 
energy laminate strength failure criterion, the minimum gage, and the 
minimum stiffness failure criterion are considered. 




ORIENTATION OF COMPOSITE PLY ANGLES 


Options and intaoT control vriabl«« .- A number of integer variables 
are included to control the flow of the program, ihese are split into 
two types: (1) those that determine the number or highest order coeffi- 

cient of various inputs; and (2) those that provide optional control for 
flow of the problem. *Ihe latter are also used to a'^oid unnecessary 
computations . 

The following are definitions of a nuniber of inputs which are under 
the control of the user: NX, NY - nund^ers and highest order of chordwise 

end spanwise polynomials, respectively. (Default: NX>3, NY*5) their 

product determines the number of displacement 2 md rotational degrees of 
freedom in the structural analysis. For the function describing the 
displacements, NX varies from 1 to NX; while NY varies from 3 to NY, 
reflecting the clamped-fixed boundary conditions at the root in the 
spanwise direction. 

NX NY 
i-l j-3 

The resulting degrees of freedom for displacement (NW) , and rotation about 

the X- and y-axes (NB ana NA) become: 

NW - NX* (NY-2) 

NA - (NX-1) *NY 

NB - NX*(NY-1) 

Description of Sample Input 

This secticn describes the output from a sample input case, it is 
intended to act as a basic guide to many of the WADES input default con- 
ditions. The actual input data and the program copy of the input data with 
all the defaults included are shown. 

Figure 23 is a direct copy of the input to the WADES program required 
to execute a minimum-weight strength design in the stand-alone program 
mode. Figure 24 is the output of the input data in figure 23. in most 
cases default values for program options have been used to demonstrate 
the minimum input required to execute the program. The input bloc)c 
designation as given in the input user guide has been written in the right 
margin at the start of each bloc)c. 


. .Half. /cni'- aV/i . iliiiri ' ‘ I 
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Block 1 is a sununary of program options. The multiple values of 
IFT, JFT,etc. , are reserved for prescribing a different number of thickness 
function coefficients when designing multilayered composite wings. Zero 
values of a given parameter usually indicate a non-active option. Block 2 
contains a summary of the variables used to define the initial geometry. 
The values for the thickness function, PTT, used in this sample data have 
been a good first estimate to the thickness distribution on a variety of 
aircraft. The value of zero for the depth function will later be replaced 
by a value computed from the alternate depth representation indicated. 

Block 3 and Block 4 contain the required material properties and 
design flight conditions. It is noted that the wing loading shown in the 
input is not used in the program, but is later computed according to the 
option, lANAL(2,j). Block 5 contains the concentrated load information 
that varies with flight condition. The value of the load at each station 
corresponds to the location of an external store or force on the wing. 

The only Block 6 constraint used is the maximum allowable stress. The 
other values are representative constraints but are not used in a strength 
design. Block 7 defines the optional parameters for the optimization 
routine, CONMIN. A zero value in this input returns to the prescribed 
default value on execution. No Block 8 or Block 9 information was read. 
These blocks are read only for the design of multilayered composite skins. 

Description of Output 

This section describes the output for the wing design progreun, WADBS. 
The particular output obtained by the user is a function of the value of 
the print control options, IPRNT and JDUMP in namelist OPTNS. The amount 
varies from nothing for IPRNT^O and JDUMP"0 to a debug level of print 

including a dximp of the various program matrices for IPRNiyo and JDUMP-5. 

A brief outline of available output including a sample case is given here. 

The basic arrangement of the output follows the calculation of the 
design. For a given 2 uialysis of the wing, and depending on the value of 
the print controls, the available output includes a printout of the 
matrices K, A, B, MX, X^MX X, and the vectors {u{, {eig{, and {x| 

for each of the flight conditions in which it is used, and a dump of 

weights, stresses, and design variables contained in the common blocks. 

The output of variables contained in common blocks is performed by routine 
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WOUT. A sample case from the F-5 design is shown in figure 23. For a 
given design cycle, the above output is optionally available for every re- 
analysis or just for the first and last designs. See the user's guide for 
the input of appropriate values of IPRNT. At the completion of the design 
cycle a summary output of common variables, and plots of the depth and 
thickness distributions, the statically loaded deformed shape and the 
normal modal shapes, the value of the flutter determinant versus Mach 
number and frequency, and the shear and moment distributions \'ersus span 
are optionally available. 

Figure 25 contains the basic summary of the current design. The 
print control required to obtain the appropriate output has been designated 
on the right side of the printout. The first line contains the title of 
the prescribed run. The root chord, semispan, and sweep emgles are a 
repeat of the input. The total wing area is the theoretical planform area. 

The weights shown contain the total weight computed from the skin and 
core weights for the designated number of wings. The skin and core weights 
are those computed by integrating the material and core distributions over 
the total wing area. The fuel-available weight is the product of the 
volume contained within the structural wing box, the density of fuel, and 
the fraction of structural volume available for fuel. The skin structure 
and core structure weights are those computed by integrating the skin and 
core material distributions over the structural planform. These structural 
weights are the weight components from which the regression analysis deter- 
mined the wing non-optimum weight factors. The locations of the centers 
of mass of the various weight components of a single wing are included 
with respect to a reference coordinate system located at the junction of 
the wing leading edge and the root chord. 

The thickness, depth, and camber coefficients are the values of the 
coefficients in the polynomial functions describing the corresponding wing 
properties. See the definition of each function to determine the power of 
^ and rj to which the coefficient is attached. If the value of the wing 
thickness-to-chord ratio is prescribed, the functional distribution for 
the depth is obtained from a linear fit of the maocimum chord depth. 

The stress distribution with respect to ^ and 11 as computed from 
von Mises ' stress resultant is printed next. The values for each of the 
design flight conditions are shown. These stresses and their locations 
are the values used to design the thickness coefficients for the isotropic 
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material wing. The component edge loadings shown next are the component 
loads per unit inch at each of the stations at which the von Mises' stresses 
are computed. These are used to design the thickness coefficients for the 
multi-layered composite material wing. The values of thickness or ply 
thicknessei as determined from the functional distribution are also given 
at each station. 

A summary of the number of computed aerodynamic and structural responses 
versus the flight conditions for which they were analyzed follows. For the 
static analysis using a constant-pressure loading the value of the angle 
of attack is only estimated from an approximate value of the lift curve 
slope. The tip deflections, natural frequencies, and flutter Mach nuinber 
and frequency (not shown in figure 25) are the values use.1 to compute the 
appropriate design constraints. 

The last item is a summary of the CPU time used during the computation 
of each of the various functions. 

Program and Subprogram Descriptions 

The following are brief descriptions of each of the routines in the 
WADES progreun. These descriptions are not intended to provide a detailed 
breedcdown of program flow, but should only indicate to the reader the 
basic usage of each routine. 

Routine Description and Comments 

Main Stand-Alone Program 

WADES Main calling program to organize reading of 

execution of analysis or design routines, 
to the appropriate routines is governed by 
parameter, ICONTR. 

Basic Structural Analysis Routines 

ANLYS Routine to organize the various analyses of the equivalent 

plate model of a trapezoidal wing. Sets up following 
solutions: 

Static: 

[Kj^+A]jwf - |Q{ , also {o|, Wj^g, w,j,g 


input and 
Branching 
choice of control 
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WINIT 


STIFF 
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Natural Frequency: 

(tK) - oi®(M)) jw( - 0 

Flutter: 

|-o)*[M] + + [Kr+AJ| - 0 

Divergence: (not availzible in ACSYNT version) 

([Kr] + X(A]){w} - 0 

Print controls IPRNT, JDUMP and control matrix lANAL are 
used to control execution of each phase of the calculation 
of the various analyses and o^:^tput. See figure 26 for flow 
chart of basic subroutine flow through routine. 

Routine to generate established boundary conditions, initial 
geometric variables, composite properties, and weights. 

Also computes approximation to depth function, d(^,r]), as 
a function of t/c and the location of meucimum t/c at 
root and tip. 

Routine to build the trapezoidal wing stiffness matrix. 
Constitutive Relation: 


r^x ^ 
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Strain-Displacement Relationship: 



[<t]|u| 


Combined Skin and Cora Orthotropic Propart ias: 


STDC3 


LOADS 


Dll - 
D22 - 
D33 - 
012 - 
D44 - 
D55 - 

Stiffness Matrix 
Freedom: 

(Kl - 


For multilayered 
through depth. 
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in terms of the Component Degrees of 
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SYM 
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composite, material properties are summed 


Routine to generate the thickness, depth and ceunber functions 
from the coefficients of polynomials. 

NTL 

tt - + t 

6-1 J J 

d ■ Wp(^,Tl) .fd(^,Tl) 

«= ■ Wp(^.n) 

Wp - (^ - TANl q) (-^ + TAN2 + RR) 

Routine to compute the work equivalent loads on a trapezoidal 
wing planform using a constant-pressure loading. Also 
solves for static equilibrium displacement vector from 
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(KrJM - |Qr( + {Q„| (Qcro} 


DPLCTN 


STRESS 


DAERO 


MASSMX 


with •xibstltution to obtain {u). 

Routine to computa the dleplacament of the leading and 
trailing edges of the wing tip. 

”le ‘ 


“tE ■ |V5TE'l-“>ri''l 

Routine to compute the von Mlses* stress and component 
strains as a function of fr«n the displacement 

vector, (u). Component skin loads at each station are 
also computed for use with multilayered composites. 

(«llt - 

(“Ik - >°'l'lk 
<»>k - ‘tl“lk 

Routine to generate the steady and unsteady aerodyneunic 
matrices A and B from piston theory. 

T 

[A] - 2tM^p^ J(l + ^ af )l‘*’wl 

Sa 

{fll - J(l + dy 

“ Sr 

Routine to generate a consistent mass matrix for the trans- 
verse inertia of a trapezoidal plate. 

/ \ 

I «1 - 5 / l^w) ( / P 

Sa ' ' 
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NA1TR2 Routlns to sot up and solve for the natural frequencies and 

eigenvectors of the free vibration of a trapezoidal plate. 

([KrJ - a)*(Ml)|w| - 0 

PLUTER Routine to compute the flutter frequency and Mach number by 

minimization of the flutter determinant. This routine 
estimates the gradients of the determinant by finite diff- 
erence and uses a Fletcher-Reeves conjugate direction algo- 
rithm to solve the minimization problem. The one-dimensional 
search is solved by ODM. 

ODM Routine adapted from Miura (ref. 16) to perform the one- 

dimensional search for the minimization of the flutter 
determinant using the golden section technique. 

FLTMIX Routine to generate the complex flutter matrix used to 

compute the determinant of the flutter equations. 

Basic Structural Design Routines 

MWT Routine to organize the analyses and the generation of 

information for the minimum weight design of isotropic 
cover sheets on a trapezoidal wing. Routine calls for 
analyses of the wing, computes gradients of the objective 
function emd initial information for geometric and 
strength constraints, and controls printing of output. 

Two versions of this routine are available: 

(1) uses the feasible direction search - CONMIN 

(2) uses a linear programming solution - SIMPLEX 

MWT43 Routine generates the gradient and constraint information 

calculations for MWT. 

FMWT Routine to generate the objective function, minimum gage, 

and strength constraints for the CONMIN version of MWT. 

FMWT2 Routine to calculate minimum gage constraints on the func- 

tional distribution of the thickness in the form to be 
used by CC^IMIM. It contains logic to by-pass calculation 
of non-active constraints. 
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MWT2 


MWT4 


MWTSTF 


PWTSTF 


FMWT3 

MWTC 


MMTC43 

Mvrrc3 

MWTC 5 
FMWTC 


Routine to compute the gradients of the objective function 
(weight of the skin) with respect to the thickness design 
variables. 

Routine to compute the location of minimum gage or strength 
constraints on a trapezoidal wing planfom. NSIG 
locations are generated using even chordwise and spanwise 
increments . 

Routine to initialize C(%IMIN parameters for the minimum 
weight stiffness and strength constrained design of isotropic 
face-sheet wings. Routine prepares control options, print 
controls, and medees initial call to COielMIN. 

Routine to organize analyses for the calculation of the 
objective function (wing weight) and the minimum gage, 
strength, deflection, frequency, flutter, and divergence 
constraints. It serves as the siibprogram called by CONMIN 
for the isotropic stiffness design, MWTSTF. 

Routine to compute stiffness (deflection, natural frequency, 
flutter, and divergence) constraints. 

Routine to initialize and organize analyses for the miaimum 
weight design of multilayered composite cover sheets on a 
trapezoidal wing. It calls for analyses, gradients of 
objective function, initial information for geomstric 
minimum gage and strength constraints, and the generation 
of output. 

Routine generates the gradient and constraint information 
for MWTC. 

Routine to initialize CONMIN parameters for the number of 
design variables and the number and type of constraints. 

Routine to initialize the minimvun gage, TL(i,2) and TLMIN. 

Routine to serve as the subprogreun for the evaluation of the 
objective function and strength, stiffness, and minimum 
gage constraints for the design of multilayered composite 
wing cover sheets. The routine is organized to compute 
the analytic gradients of strength and min.imum gage con- 
straints. It is also set up for the calculation of finite- 
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mnc3 


MWTCST 


PWTCST 


SWITCH 


SWTCHS 


SWTCHC 


SWTCHD 


diff*r«nc« gradisnts of displacoinont, froquonq^, and 
fluttar constrainta. This la an adaptation of routina 
COMPOS. 

Routina to computa tha analytic gradlanta of tha gaomatrlc 
minimum gaga constraints with raspact to tha thlcknass 
function daslgn varlablaa and store tham In tha A matrix 
In common block CNMM2. 

Routina to Initial Iza COMMIN paramatars for tha minimum 
walght daslgn of multllayarad composlta covar shaats to 
satisfy stlffnass and strangth constraints. It proparas 
tha control options, print controls, and initializations, 
calls CONMIN, and controls tha ganaratlon of output. 

Routina to sorve as the subprogram for tha calculation of 
I'ha objective function (wing walght) and strangth, minimum 
gage, deflection, frequency, flutter, and divergence con- 
straints for the design of multilayered composite cover 
sheets of a trapezoidal wing. Routine calls FMWTC to 
calculate constraints. 

Routine to store and retrieve analysis control parameter, 
lANAL. lANAL is stored and retrieved as a function of 
flight condition when calculating the analyses corresponding 
to active constraints. 

Routine to test for active Isotropic strength constraints 
and to set appropriate control parameter, lAHAL, on or off. 
This Is used only during stiffness design by FWTSTF. 

Routine to test for active composite strength constraints 
and to set appropriate analysis control parameter, lANAL, 
on or off. This is used only during stiffness design by 
FWTSTC. 

Routine to test for active stiffness (deflection, frequency, 
flutter or divergence) constraints and to set appropriate 
analysis control pareuneter, lAHAL, on or off. 
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Input/Output Routinaa 


WIN PUT 

WOUTPT 

WOUT 

MATRLS 

C0MP14 

COMP25 

COMP30 

MXOUT 

DISTRB 

PLTDPL 

PLTPLT 

PLTLDS 


Routine to read and make a copy of the input. The first 
pass through this routine sets the program defaults. See 
WADES user's manual for a description of variables read. 

Routine to organize printing of detailed output and plots. 
Amount of output is determined by choice of control param- 
eters IPRNT, JDUMP, and I PLOT. 

Routine to output all variables in common blocks BWSAV 
and BSTRCW. 

Routine to initialize material properties for three different 
wing cover sheets: (1) aluminum; (2) titanium; and 

(3) graphite/epoxy. 

Routine to read the analysis/design variable transformation 
for composite material thickness variables of ply orienta- 
tions. 

Routine to read and write a copy of the composite input 
data for direct input of Izunina properties. 

Routine to read and initialize the additional optional 
parauneters of CONMIN. Program is initialized to default 
values on first call. 

Routine to generate the output of a general matrix, A(N,M). 

Routine to produce a contour plot of the thickness, depth, 
or caunber functional distributions versus (x,y) on the wing 
planform. 

Routine to produce a contour plot of the lateral displace- 
ments or mode shapes versus (x,y) on the wing planform. 

Routine to produce a contour plot of the flutter determinant 
versus u> and Mach number. 

Routine to generate plots of the loading per unit span due 
to the weight of the skin and core and the aerodynaunic 
spanwise loading. An output of the area and moment of 
inertia of the structure is also included. 


PLT8HR Routln* to gonorato • plot vorsus upon of tho shoar distri- 

butions dus to skin, cors, and aarodynamic loadings on a 
trapazoidal wing. Tlta rasulting banding-roomant distributions, 
including discrata loads, arc also output. 

C0MP19 Routine to output a table of active constraints during 

composite design. 

General utility Routines 

l^ese are general routines for matrix analysis and algebraic function 
manipulation. Those routines which operate on functions of (x,y) assume 
the polynomial function is of the form 

P(x,y) " «nd j-l,JP 

where i-fj-2 < max(IP,.7P) -1 , and where the coefficients P^^ are stored 
as a single vector. For example, if ZPb 3 and JP"4, the polynomial is 
written 

p(x,y) ■ Pj + P^x + P^x® + P^x® 

+ PgY + PgXy + PgX^y 

Description and Comments 
Routines to multiply two polynomials together: 

C(^,Tj) - A(4,ri)x 

ADOC, Routines to add two polynomials together: 

AODC2 

C(?,ri) » A(^,ti) + 

OXZC Routine to generate the derivative of a polynomial with 

respect to ^ : 

C(4.Tj) ■ dA(^,Tl)/d^ 

VALUE Routine to evaluate the polynomial fer a given value of : 

RESULT - P(^)^) 


Routine 

MULTC, 

MULTC2 
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Function to evaluate a polynomial of two variables for a 
given value of • 

VALUE2 - 

Routine to evaluate a polynomial of two variables along a 
line of comtant to form a polynomial function only 

of ^ ; 

AX(^) - 

Routine to generate a set of polynomials, their first two 
derivatives, and the product of any two ccr'^inations of the 
above. These are used as the chordwise displacement func- 
tions. Both ordinary and orthogonal polynomials are 
available. 

Routine to generate a set of spanwise polynomials, their 
first two derivatives, and the product of any two combi- 
nations of the above. These are used as the spanwise 
displacement functions. Both ordinary and orthogonal 
polynomials are availedjle. 

Routine to perform the double integration over the wing 
surface of the product of three polynomials: 


RESULT - R*S j H(tj)D(^,t 1) P(^)d^ dt) 


where H and P are products of the displacement functions 
and D is a function distributed over the wing surface. 

Routine to generate a table of integrated values of 
T|^”^ over a given trapezoidal planform: 


BINTG(i,j) - R»S 


if 


d^ dTl 


Function to evaluate the planform polynomial: 
FPLAN - - TANl (-?j^ + TAN2 t\^ + RR) 
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PXI 


SOLVBl 


EI6ENR, 

TRED2, 

TQL2 


PAS003 


PAS030 


WDUMP 


CPUTI.M 


BUFFER 


Routine to compute the ^ location in the nondimensional 
Cartesian coordinate space of a constant fraction of chord 
on the trapezoidal wing planform: 

FXI - TANl + C(RR - (TANl - TAN2) 1 

Routine to solve a system of equations by Gauss elimination; 

A.X - b 

Routines to solve the real eigenvalue problem in the form: 
A*X - XB*X - 0 

where A is symmetric and B is symmetric positive 
definite. 

Routine adapted from program PASS to reduce a positive 
definite FxN matrix contained in A to an NRxNR matrix: 

(All - A12*A22”^A21) A12 
(A22“^A21) 1 


A - 


All A12 
A21 A22 


Routine adapted from program PASS to calculate the absolute 
value of the determinant of a complex NxN matrix by Gauss 
elimination without pivot search. 

Routine to putput the complete set of common blocks of the 
WADES progrzun *n the event of a major error return from a 
subprogram. Progrzun is terminated in case of error. Routine 
may also be ussd to obtain a dump of the common blocks in 
certain cases. ACSYNT version of routine only sets the 
return error code. 

Routine to compute the 2d3Solute and relative CPUTIM during 
execution. 

Connects the analysis and design variables in routines MWT 
and MWTSTF. 


Composite Analysis and Design Routines 

These routines are obtained directly or adapted from the COMPOS 
progreun. 
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COMP07 

COMPOS 

C0MP09 
COMPIO 
COMP 11 
COMP 12 
COMP17 
COMP18 
COMP 20 
COMP21 
COMP26 
COMP27 
COMP28 
COKP29 


Computes the transformation: 

[Oj] 

Sums the composite stiffness from Izunina properties: 

NL 

W - E 

computes gradients of composite strain with respect to 
Icunina thicknesses; 

Computes gradients of strength failure criteria with respect 
to lamina thicknesses; 3|g}/dt^<. 

Transforms ply or thickness design variables into analysis 
variables by appropriate transformation. 

Transforms thickness design variable of polynomial function 
into respective leunina thicknesses. 

Routine to multiply gradients by analysis/design variable 
transformation and store in A matrix of CCWMIN. 

Multiplies gradients by analysis variable/ thickness function 
transformation and stores in A matrix of CONMIN. 


Computes the gradients of the composite strain with respect 
to the l€unina ply orientations; 3(€g}/dt9. 

Computes gradients of strength failure criteria with respect 
to the lamina ply orientations; d|g)/d(9. 

Computes constituent properties and makes a summary output 
of composite and leunina properties. 


Computes gradients of lamina transformations with respect 
to ply orientations; 

Computes gradients of composite stiffness Aj, with respect 
to thickness variables; dAc/dt^. 


Computes gradients of composite stiffness with respect 

to ply orientation; dAf./dO, 


Computes composite failure criteria, |g}. 


COMP 31 


Character Plotting Routines 


ARPLOT 

PL0TA2 

PL0TA5 

PL0TA6 

PL0TA7 

PL0TA8 

PL0TA9 


Routine generates a three-dimensional character plot of 
vectors of X, Y, and Z. 

Routine to plot NY curves of Y versus X. 

Routine to provide rounding of scaling varii^les to a 
specified exponent. 

Routine to round off scaling to nearest acceptable plotting 
scale. 

Routine to initialize a row of plotted output including 
location of X and Y axes. 

Routine to select scales, round maximum and minimum values 
to acceptable values, and locate X and Y axes. 

Routine to generate linear and logarithmic scaling of the 
Z coordinate variables to accept 2 d»le values for use with 
contour plotting. 


Numerical Optimization Routines 


CONMIN 


SMPLXl , 
SMPLX2 


Subprogram and associated routines for the solution of 
constrained minimization problems. Called by MWT, MWTSTF, 
MWTC, and MWTCST. See reference 14 for user's manual and 
description of parameters. 

Subprograms for the solution of the linear progrcunming 
problem. Routine is called by alternate form of MWT. 


Use of Modified Strip Loading Routines 

In order to use the codes for the modified strip loadings described 
in Appendix C ^wo routines used to generate loads in the current WADES 
program have to be replaced with their corresponding equivalents before 
the start of execution. The codes for the strip loads are arranged such 
that no program changes have to be made in order to accommodate them in 
the present program flow. The two routines to be substituted for are 
DAERO and LQADi^ . Seven additional routines are also called by the above 
two substitute routines. 
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Since the modified strip analysis is limited to subsonic flight 
conditions, the use of the corresponding computer codes also restricts the 
progreun to the consideration of Mach numbers perpendicular to the quarter- 
chord line of less than unity. That is, the free-stream Mach ntimber must 
be such that 

cos A ^ ^ 

Though this is an adjsolute limit on the program, practical considerations 
suggest that Mq cos A should not be in the transonic range either. This 
Mach number restriction has also led to the elimination of unsteady aero- 
dynamic codes using piston theory, and thus the calculation of supersonic 
flutter is no longer a progrsun option in this mode. The only input modi- 
fication required to use the modified strip analysis code is for the 
generation of the spanwise shear and moment plots. The plotting control 
option must be changed to IPL0T(8)>3 to obtain plots. All associated 
panel geometry and transformations required to interface the discrete 
loadings with the plate model are generated internally. 

Description of Modified Strip Loading Routines 

The following are brief descriptions of each of the modified strip 
analysis routines in the WADES program. These routines are inserted at 
progrzun load time. 

Rout ine Description and Comments 

DAERO Alternate routine to compute and invert aerodynamic matrix, 

[Sj^], and to generate the transformation (T^l and [T^l 
for modified atrip static load analysis. Transformation 
[T^] defines the relations between the local amgles of 
attack, !aj, and the lateral displacement, (w}. Transfor- 
mation fTwl defines the relation between the work 
equivalent load vector and the strip panel lift. The 
matrices are stored as: 

A(i,j,l) -['4m^J[Sj"S A(i,j,2) - [T^] , and 

A(i,j,3) - [Twl 
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LOADS 

EQUILK 

MXSl 

PLOAD3 

SHEAR 

PASOOl 

PAS002 

PASINV 


Substitute routine to compute the work equivalent modified 
strip loads using a Ritz-type analysis and to solve for the 
static displacements vector, (u). Spanwise shear and moment 
distribution plots are called for here because of data 
transfer limitations. 

Routine to compute the values of the root angle of attack, 
tail load, lift distribution, trim aileron deflection, and 
structural deformation required for static equilibtium in 
maneuvering flight for symmetric and nonsymmetric planforms. 
Routine has been adapted to accept the stiffness form of 
structural and aerodyneunic matrices. 

Routine to compute the aerodynamic influence coefficient 
matrix, SI. This matrix is computed for symmetric and 
nonsymmetric planforms and for symmetric and antisymmetric 
loading conditions. For symmetric planforms, y>0 is the 
centerline, and the plane of symmetry is the X-Z plane. 

Routine to generate plots of the loading per unit span due 
to the weight of the skin and core and due to spanwise 
loading from modified strip theory. An output of the area 
and moment of inertia of the structure is also included. 

Routine to generate spanwise plots of the shear distributions 
due to skin, core, and modified strip aerodynamic loadings 
on a trapezoidal wing. The resulting bending -moment dis- 
tributions, including discrete loads, are also output. 

Routine adapted from the PASS program to form the LU 
decomposition of the positive definite matrix A. 

Routine adapted from the PASS program to solve the system 
of equations LU*X •> B, by forward and backward substitution. 

Routine adapted from PAS002 to form the inverse of matrix A 
by forward and backward substitution on the identity matrix. 
The forward substitution initializes the identity matrix. 
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Use of Supersonic Piston Theory Loads 

In order to use the codes for the supersonic piston theory loadings, 
only the LOADS routine in the current WADES program has to be replaced 
with its corresponding equivalent. This direct substitution is made prior 
to start of execution. The codes for the piston theory leads are arranged 
such that no program changes have to be made in order to make the substi- 
tution. This LOADS routine uses the information from the default routine 
DAERO directly. 

The theory in this routine is generally limited to aerodynamic flow 
fields in which supersonic Mach numbers exist everywhere on the wing 
surface. The aerodyneuk;ic pressure distribution is defined in the section 
on loads in the WADES user's manual. The equations solved to obtain the 
deformed shape of the wing are: 

Vertical equilibrium: GLR - Ls^(ao) - 

Matrix equations: [Kr+a]|w| - (Os^^a^)} +|0w)+ lOcml 

The two equations are solved iteratively to obtain {w} and a^. The 
remaining displacements, {u}, are determined in terms of {w}. 

Implementation on the Computer 

The WADES program has been written in standard ANSI FORTRAN IV 

language. It has been run on both the COC 7600 and IBM System 360/67 

computers. Though the language is standard FORTRAN IV, minor changes in 

the codes must be made in order to make the program operational on both 

systems. These generally have to be made to account for the machine- 
dependent functions of the two computer systems. A summary of the diff- 
erences between the codes for the two computers follows. 

(1) Because the calls to the system clock vary between computers, 
the routine CPUTIM must be modified to call the machine-dependent clock 
of the particular computer. 

(2) Routine MXOUT passes an eight-character title of the matrix to 
be printed. Because cf the differences in word lengths between the two 
computer system, MXOUT must be modified to accept either two four-byte 
words for the IBM 360 or one ten-byte word for CDC by changing the 
dimension of FAR. 
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(3) Bttcause of the differences in computer word length and the 
resulting loss of accuracy for the smaller word size, it is recommended 
that routines performing the polynomial integration be converted to double 
precision. All lines requiring switch changes for precision in the presei.t 
version contain CDC or CDBL in the first four columns. Removal of the 
comment character on the appropriate double precision statements indicated 
is sufficient. 

(4) Alternate routines for use with the CDC overlay systems were 
written in a number of cases. The non-overlay routines should be used 
with the IBM 360. 


CONCLUSIONS AND RECOMMENDATZC»IS 

The primary objectives of this study were the incorporation of the 
WADES programs as a wing design module of the ACSYNT progreun and the cor- 
relation of the weights and material distributions with existing aircraft. 
The first objective, integration of the WADES progreun into ACSYNT, was 
accomplished, and a demonstration case using the F-5A was executed to 
demonstrate progreun convergence. Default values for a nundaer of optional 
pareuneters were defined and programmed to increase the progreun useedsility. 
Similarly, interfaces were written for variables used to define the wing 
geometry and loads in terms of ACSYNT descriptors. 

The second objective, the correlation of weight and material distri- 
butions with existing aircraft, was investigated with two simultemeous 
approaches. In the first, a detailed comparison of material distribution 
and loading conditions was undertedcen. This included a surface fit and 
detailed analysis of the thickness distribution of the upper skin of the 
F'*5A. Comparisons of the assumed design loading conditions with the actual 
values were shown for maximum symmetric pull-up, and landing. The data 
were then recomputed using the design algorithms defined in the WADES 
progreun. In the second approach, the statistical determination of a basic 
non-optimum weight factor was undertedcen using data from U. S. Air Force 
and U. 3. Navy fighter aircraft. This included the analysis and design by 
the WADES program of the aircraft wings. A regression analysis using the 
actual and estimated values was made to correlate results for several wing- 
weight approximations. 
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It is concluded that the comparisons demonstrated in the F-3 study 
show reasonably good correlation for the estimation of the loads and material 
distributions, but do not show strong statistical correlation for the break- 
down of weights considered. The subsonic loads in figure 10 show the best 
correlation with actual F-5 data, especially for torsion. The results of 
the surface fit of the thickness distribution show good correlation with 
the deformed wing shape in figure 13, but they indicate a program deficiency 
in modeling the structural planform. Though the function fits the thick- 
ness distribution closely, the inability of the structural model to repre- 
sent the discontinuity in the stiuctural planform results in grossly 
underestimating the material distribution at the tip (fig. 5) . Conversely, 
the lack of modeling of the wing-fuselage junction produces the excessive 
bending moment at the root and a major discrepancy in the loads in fig- 
ures 9 and 10. The resulting design of the F-5 wing by the WADES program 
similarly overestimates the required material distribution at the root 
as indicated in figure 5. 

Though the material distributions do show some correlation with actual 
F-5 values, the calculated weight, even when the functional fit of the 
data is used, are typically 25 to 50 percent low on the weight of the 
structure. The even higher coefficients obtained in the regression 
analysis in Tedsles VI and VII, which also reflect the additional weight 
of leading-and trailing-edge structure, indicate the relatively high pro- 
portion of non-optimum structural weight. The values of the non-optimum 
weight factors recommended for usage with the present program are found 
in the sections on regression results for isotropic and composite wing 
designs. 

As a result of the present correlation activity, a number of defi- 
ciencies and omissions were identified. They include deficiencies both 
in geometric and structural modeling and in critical capabilities for 
minimum acceptcdsle aeroelastic modeling. The recommended corrections and 
improvements that can be made within the framework of the current progr2un 
are as follows. 

(1) The current modeling of the wing-fuselage junction is inadequate 
as seen from the loads in figures 9 to 11. The assumption that the wing 
is clamped at the centerline rather than at a finite body radius can 
seriously misrepresent the spanwise shear and moment distributions, 
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•apaclally for low-.ispact-ratio configurations whsra a significant portion 
of the planform area is contained within the body. This effect is twofold: 
first, the aerodynamic load is considered to be distributed over the theo- 
retical planform rather than the exposed planform; and secondly, the 
increased structural ipan produces higher bending moments and therefore 
higher required material gages at the root. The modifications recommended 
to correct this modeling deficiency are first to use two concentrated 
loads located at the wing-fuselage junction to react the shear and torsion 
at the root and secondly to reformulate the statics to distribute the 
pressure load only over the exposed planform. 

(2) Over half of the F-5 wing panels were identified as buckling 
critical as seen from the stresses in Ted)le B-II in Appendix B. The only 
strength failure criterion currently used with the WADES program is an 
nltlmate failure criterion on the material itself. The modification 
recommended to correct this deficiency is to include a simplified buckling 
failure criterion such as the ones in figures B-1 and B-2. This could be 
effected by modification of the failuxa criteria in the design phase to 
include a buckling stress interaction curve. 

(3) The current structural model of the skin is that of an isotropic 
covei plate which ignores the effects of stiffeners and spar flanges. It 
is recommended that this be reformulated as an effective skin-thickness 
model. The effective thickness would then be computed on the basis of 
stiffener spacing and size. This modification would complement the inclu- 
sion of the previously recommended buckling constraint by establishing the 
buckling pemel size. It will also be used to improve the estimation of 
non-optimum weight. 

(4) The estimation of the cross-sectional area of structural material 
in figure 5 identified a prograum deficiency, naunely, that the restriction 
of the structural planform to a single trapezoid severely restricts the 
modeling of the actual plamform. The mismatch is most obvious in figure 2 
and is similarly reflected in the estimated structural weights. The modi- 
fication recommended to improve this problem is the division of the plan- 
form into three trapezoidal segments. This should allow the flexibility 
to include a wing carry-through structure and the cutout for the landing- 
gear structure. It may be implemented by modification ol the integration 
tables and the geometry necessary to describe them. 
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(5) Th« correlation coefficients determined from the regression 
analysis for the group of U. S. fighters considered did not show a con- 
clusive trend for the combinations of variables considered. A further 
breakdown of the weights computed by the WADES program should be made. It 
is recommended that the regression on the weights be repeated with more 
attention toward obtaining correlations of component weights as well as 

of the total weight. 

(6) The integration of chs WADES program required the generation of 
the critical structural design loads from the flight profile of the aircraft 
mission. The temporary implementation of this was to design a specific 
routine to compute the critical conditions baaed on previously determined 
requirements of the aircraft, it is recommended that this routine be more 
generalized and that a check of FAR-25 or appropriate MIL requirements be 
made to determine a critical flight profile. 

(7) The preliminary flutter analysis indicated that supersonic piston 
theory availedsle in the current program version was not appliced}le for the 
flight profile of the F-5. It is recemunended that an availedsle subsonic 
flutter calculation be added to the program to fill this gap. 


NIELSEN ENGINEERING & RESEARCH, INC. 
Mormtain View, California 
August 1975 


TABLE I 


SUMMAFY OP P-5A/B WING GEOMETRY 


Aerodynamic Planform Oeacription 


Engineering 

Symbol 

Fortran 

Name 

Value 

(t/C) 

TCR 

0.048 

(V') tip 

TCT 

0.048 

x/c ( t/c ) 

XTCR 

0.50 

°root 

R 

134.5 in 

b/2 

SPAN 

151.5 in 

®LE 

THETl 

32.0° 

9te 

THET2 

-5.0° 

Structural 

Planform Description 

Engineering 

Symbol 

Fortran 

Name 

Value 

(x/c)le 

XLE 

0,15 

( 1 ” X/ C ) 

XTE 

0.55 

PVA 

PVA 

0. 



TABLE II 


SUMMARY OP THREE CRITICAL FLIGHT COMDITIONS 


Northrop Identification; 

#104 

#356B 

#123C-5 

Maneuver Description*; 

Symmetric 

Pull-Up 

Dynamic 

Landing 

Symmetric 

Pull-Up 

Engineering 
Symbol (\mits) 

Fortran 

Name 

Flight Conditions 

1 

2 

3 

h (ft) 

ALT 

21,500 

0 

0 

Mo 

XMO 

1.05 

0.22 

0.90 

SGW (lbs) 

— 

11,543 

12,200 

11,591 

Wfuel (lbs) 

WPUEL 

0 

0 

0 

'^wing (lbs) 

WTWING 

1041.7 

1041.7 

1041.7 

NZ^lt 

ANZ 

9.8 

4.0 

9.8 

Mbody (lbs) 

WBODY 

9,597 

10,432 

9,359 

XcG (% MAC) 

PMAC 

11.4 

10.2 

10.7 


Concentrated Mass Weight Information (lbs) 


Description 


Flight Conditions 


Location 


x/c 


Landing Gear (main) 
AIM-9B/Tip Tank 
Pylon WS 85 
Pylon WS 114 


248 

204 

0 

0 


-4327** 

115 

0 

0 


248 

115 

112.75 

119. 


0.55 

0.10 

0.30 

0.30 


0.40 

1.00 

0.;j6 

0.75 


** 


Speed brakes are closed during all flight conditions. 
Percent of lift ac landing equals 100 percent. 
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TABLE VI 

REGRESSION ANALYSIS ON WEICIlt PREDICTION 01 <!A 'ES PROGRAM 


CCNSTANT M/S LOADS ' OBJ-MTSS 


AIKCRAFT 

MWING 

MTSS 

ntcs 

S 

F-5A/B 

1041.7 

381 .1 

221.8 

24403.0 

F-5A/0 

1041. 7 

318.8 

221.8 

24403.0 

A-6 A 

4778.2 

1475.8 

3318.3 

73766.0 

A-7A 

3275.8 

1145.9 

1736.8 

53944.0 

F-<rC 

46 71.0 

1682./ 

1796.9 

76607.0 

A4D-2N 

1587.3 

322 .5 

1095. 1 

37293.0 

F8U-1 

2765.8 

965.1 

1352.6 

53147.0 

F'yF-b 

2600 . 1 

679.6 

1554.9 

42814.0 

F-104A 

1 180.6 

340.3 

409.8 

30612.0 

F-1 C5B 

3404. o 

2843 .8 

879.6 

53310.0 


SUMMARY OF REGRESSION RESULTS 


CONSTANT m/s LOADS - OBJawTSS 

WlaB( I l*MTSS*B(2 )*WTCS 

w2aB2»WTSS<-MTCS 

M3»B3*MTSS 

h‘»=BS(l)aWTSS*BSI2MS/288. 


B( 11 = 
PS (1) a 

3.9971 

0.3269 

B( 2) a 
B$(21a 

1.1603 

14.6781 

B2a 

1.1680 

83 a 

2.1159 


AI BCRAl^ T 

WMlNG 

Ml 

PI 

M2 

P2 

m3 

P3 

M4 

F-5A/L' 

1041 .7 

641.8 

>38.4 

666.9 

-36.0 

806.4 

-22.6 

1368.3 

F-5A/B 

1041.7 

579.7 

-44. 4 

594. 1 

-43.0 

674.5 

-35.2 

1347.9 

A-6A 

4 778.2 

5388.0 

12. 8 

5042.0 

5.5 

3122.6 

-34.6 

4241.9 

A-7A 

3275.3 

3192.5 

-2.5 

3075.2 

-6.1 

2424.6 

-26.0 

3123.9 

F-4C 

4fe71.0 

3798.6 

-18. 7 

3762.2 

-19.5 

3560.4 

-23.8 

4454.4 

A40-2N 

1587.3 

1614. 1 

1.7 

1471.8 

-7.3 

682.4 

-57.0 

2006.1 

F8U-1 

2765.8 

2558.7 

-7.5 

2479.8 

-10.3 

2042.0 

-26.2 

3024. 1 

F9F-6 

2600. 1 

2512.8 

-3. 4 

2348.6 

-9.7 

1437. 9 

-44.7 

2404.2 

F-104A 

1180.3 

823.0 

-30.3 

807.3 

-31.6 

720. C 

-39.0 

1671 .4 

F-1C5B 

3404.0 

3873.6 

13.8 

4201.0 

23.4 

6017.1 

76.8 

3646.6 
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TABLE VII 

REGRESSION ANALYSIS ON WEIGHT PREDICTION OF WADES PROGRAM 


MUOIFICD STKIP THEORY LOADS ( TN3030I - OBJ«MrSS 


AIRCRAFT 

WwING 

WTSS 

WTCS 

S 

A-6 A 

4778. 2 

942.2 

3318.3 

73 766.0 

A-7 A 

3275.8 

863 .6 

1736.6 

53944.0 

F-4C 

4671 .0 

1156.4 

1796.9 

76607.0 

F8U«l 

2 765.8 

778.9 

1352.6 

53147.0 

F9F-6 

26Gl> . 1 

376.1 

1554.9 

428U.0 

F - 1 A 

1180.8 

223 .0 

409.8 

30612.0 



SUMMARY OF 

REGRESSION 

RESULTS 


‘'nOIFIFD STRIP THEORY LOADS 5TN3C30> - OflJ«RTSS 


Wl»B(l)»»i(TSS+8(;i)*WTCS 

w2»tt20WTSS^WTCS 

w3*B3’*WTSS 

W4 = aS( I J»WTSS«-BS (2 »*S/288. 


B( 1 ) = 

2.5891 

B( 2) » 

0. 7553 

B2= 2. 0416 

83* •♦.2795 

6S ( 1 1* 

0.7588 

BS(2)» 

14.2054 




AIRCRAFT 

MtalNG 

W 1 

PI 

w2 

P2 

83 

P3 

N4 

P4 

A -6 A 

4778.2 

4945.9 

3.5 

5241.9 

9.7 

4032. 1 

-15.6 

4353.4 

-8.9 

A-7A 

3275.8 

3547.8 

8. 3 

3499.9 

6.8 

3695.8 

12.8 

3316.1 

1 .2 

F-^C 

4671.0 

4351.3 

-6.8 

4157.8 

-11.0 

4948. 8 

5.9 

4656.1 

-C.3 

FOU-1 

2765.8 

3C38.3 

9.9 

2942.8 

6.4 

3333.3 

2C.5 

3212. 5 

16.2 

F9F-6 

2600. 1 

2148.2 

-17.4 

2322.7 

-10.7 

1609.5 

-38.1 

2397.2 

-7.8 

F-104A 

1 180.8 

886.9 

-24.9 

865. 1 

-26.7 

954.3 

-19.2 

16 •? 1 

42.2 
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TABLE VIII 

REGRESSION ANALYSIS ON WEIGHT PREDICTION OF WADES PROGRAM 


COMPOSITE 

DESIGN - 

CONSTANT 

W/S LOADS 

- 701 WMING 

IKCRAFT 

WM ING 

WTSS 

WTCS 

S 

F-5A/B 

729.0 

89.0 

221.8 

24403.0 

A-6A 

33*»5.0 

368.5 

3318. 3 

73 766.0 

A-7A 

229C.0 

382.6 

1736.8 

53944.0 

F-4C 

3270.0 

A74.8 

1 796.9 

766O7.0 

A<.D-2N 

1 IIO.O 

3A1.8 

1095.1 

37293.0 

rau-i 

1935.0 

284 .6 

1352.6 

53147.0 

F-ICAA 

825.0 

191 .6 

409.8 

30612.0 


SUMMARY OF REGRESSION RESULTS 

COMPOSITE DESIGN - CONSTANT K/S LOADS - 70* WWING 

wi3d(ii*wrss»Q(2i-wrcs 

n(2=B2«WTSS«-l*TCS 

WJ=B3»WTSS 

W^-OSm^taTSS^BS (2 )»S/288 . 


B( 1 ) = 

3.1780 

8(21 = 

0. 

6 750 

82= 1.6218 

83= 6.^188 

BSU > » 

-1 .00 15 

BS(2» = 

14. 

5329 




AI RC PAFT 

WW ING 

WI 

PI 

W2 

P2 

m3 

P3 

M4 

P4 

F-5A/B 

729.0 

432.7 

-40.6 

360.1 

-49.8 

571.3 

-21.6 

1083.5 

48.6 

A-6A 

3345.0 

3412.9 

2. 0 

3915.9 

17. 1 

2305.3 

-29.3 

3110.1 

-7.0 

A-7A 

2290 .0 

2389.1 

4.3 

2357.2 

2.9 

2455. 6 

?.2 

2386.5 

-8.9 

F-4C 

3270.0 

2722.9 

-16. 7 

2566.9 

-21.5 

3047.7 

-6.d 

3076.8 

-5.9 

A40-2N 

1 1 10.0 

1826. 1 

64. 5 

1649.4 

48.6 

2194.0 

97.7 

1314.0 

IB. 4 

F8U-1 

1935. C 

1818 .3 

-6.0 

1814.2 

-6.2 

1 826. 8 

-5.6 

22C9.0 

14.2 

F-104A 

825.0 

885.8 

7.4 

720.5 

-12.7 

1229.8 

49.1 

1226.<« 

48.7 
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TABLE X 


OPTIONAL MATERIAL PROPERTIES DEFINED IN MATRLS 
BY CHOICE OP PARAMETER, ITYPES 


ITYPE 

1 „ 

2 

Material 

2024-T3 Aluminum 

6A1-4V Titanium 

Ex 

10,500,000 psi 

16,000,000 psi 

Ey 

10,500,000 psi 

16,000,000 psi 

Gxy 

4,000,000 psi 

6,200,000 psi 

Gxz 

4,000,000 psi 

6,200,000 psi 

Gyz 

4,000,000 psi 

6,200,000 psi 

vxy 

0.3 

0.3 

vyx 

0.3 

0.3 

SmaxT 

63,000 psi 

117,000 psi 

SmaxC 

63,000 psi 

117,000 psi 

P 

0.100 Ibs/in® 

0.160 Ibs/in^ 


^Reference: MIL-HDBK-5, August 1962. 

Valid to 20''°P for 1000 hours. 


Composite Material Launina Properties 


ITYPE 

3 

Material 

Graphite/Epoxy* 

ELll 

21,000,000 psi 

EL22 

1,700,000 psi 

GL12 

650,000 psi 

GL13 

650,000 psi 

GL23 

650,000 psi 

vll2 

0.21 

vl21 

0.017 

pi 

0.056 Ibs/in^ 

SLUT 

180,000 psi 

SLllC 

180,000 psi 

SL22T 

8,000 psi 

SL22G 

30,000 psi 

SL12S 

12,000 psi 

SL23S 

12,000 psi 

PH ID EL 

0.0104 


*Single-ply thickness, ■ 0.008 in; 

void fraction, ■ 0.0; and 
fiber fraction, k^ ■ 0.50. 


















Spanwise direction » y (inches) 



Chordwise direction, x (inches) 


Figure 3.- Contour plot of thiclcness distribution, t, obtained 
from least-squares fit of F-5 upper wing cover plate. 

I 


I 



opcinwi.se axirection, y ixncnesj 


T 


Contour values are 
in inches 



Figure 4.- Contour plot of effect!^ 
obtained from least-squares fit ( 



n, X (inches) 

/e thickness distribution, te 
jf F-5 upper wing cover plate 


. : ir- rff 




Area, square inches 







Figure 6.- F-5A/B n« 3 ment of inertia, Ix, of structural material versus span 
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(AMFLT^j) 


Figure 7 




- Static balance of forces and pressures on wing 










Torsional moment (10 in- lbs) 


Condition 

123C-5 

.■j 

M - 0.90, 

h - S.L. 

1 

Mz - 6.5, 

Wt • 11,591 lbs 

j 

Empty tip 

tan)c8 and pylons only 


at WS 85 

and WS 114 


— 

• F-5A/B (Northrop, 
ref. 7) 



" WADES constant W/S 



pressure loading 

iC 

i 

t 


Wing station, inches 


> 




xgure 9.- Ultimate spanwise wing loads resolved to 35-percent 

chord axis. 







condition 123C-5 
M - 0.90, h - S.L. 

Ng - 6. 5, Mt - 11,591 lbs 
Empty tip tanks and pylons 
only at WS 65 and WS 114 


P-5A/B (Northrop, 
ref. 7) 


WADES modified atrip 




Bending moment (10 in- lbs) 







Figure 12.- F-5 wing vertical deflection for test limit load 

wing condition 104. 
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ACSYWT 



Figure 14.- WM)ES/ACSYNT subroutine floi# chart. 
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The equivalent values in terms of ACSNT program parameters are: 


Description 
Root chord 
Tip chord 
Semispan 

Reference sweep angle 
Leading-edge amgle 

Trai ling-edge angle 


WADES ACSYNT 

R ROOTWG*SPANWG 

CTIP TRWG*ROOTWG*SPANWG*SCALC(NCHORD) 

SPAN SPANWG/2. 

A A 

0 ARCTAN [TAM/^EWING (NCHORD) *CTIP/SPAN 

+ EWING ( 1) *R/S PAN) 

9 ARCTAN [TAN 0 , + ( CTIP- P ' /S PAN ) 

2 i 


Figure 20.- WADES/ACSYNT wing geometry. 














STRUCTURES MODULE 
WADES (WIND ANALYSIS AND DESIGN) 
VERSION DATED IS SEPTEMBER 1975 


THE FOLLOWING IS X BRIEF DESCRIPTION OF THE INPUT NECESSARY FOR 
OPERATION OF THE WADES STRUCTURES MODULE WITHIN ACSYNT. THE DEFAULT 
VALUES ARE INDICATED IN PARENTHESES. WADES USES A SYSTEM OF CONSISTENT 
UNITS INTERNALLY. THAT IS. ALL UNITS USED THROUGHOUT THE PROGRAM ARE 
DEFINED AND MAINTAINED AS THOSE UNITS WHICH ARE SPECIFIED AT THE TIME OF 
INPUT. BECAUSE AVAILABLE MATERIAL PROPERTY VALUES ARE GENERALLY IN 
INCHES AND POUNDS THIS SYSTEM HAS BEEN ADAPTED FOR OPERATION WITH 
ACSYNT. ACCORDINGLY. THE PROGRAM INTERFACE ROUTINE CONVERTS ALL ACSYNT 
DIMENSIONAL VALUES TO INCHES FOR INVERNAL USE. THEREFORE. ALL FOLLOWING 
REFERENCES FOR LENGTH ID AND WEIGHT 'fl WILL USE INCHES AND POUNDS 
RESPECTIVELY. 


I TITLE I ALPHANUMERIC HEADING FOR OUTPUT FORMAT I 19AA) I 


1 BLOCKl I NAMELIST COPTNS 


I 


DESCRIPTION 


I DEFAULT I 


LOPTNS 

IPRT • INPUT PRINT CONTROL OPTION ID) 

0 • 00 NOT PRINT COPY OF INPUT DATA 

1 > PRINT COPY OF INPUT 

IPRNT > OUTPUT PRINT CONTROL 10) 

0-00 NOT PRINT OUTPUT 

1 » PRINT FIRST AND LAST DESIGN OUTPUT 

2 * PRINT OUTPUT ONLY DURING 1-0 DESIGN STEP 

3 • PRINT ALL OUTPUT DURING DESIGN 

JOUMP • DEBUG LEVEL PRINT CONTROL ID) 

0 • DO NOT PRINT OUTPUT 

1 ■ PRINT MINIMAL WEIGHTS AND CURRENT DESIGN 

2 ■ PRINT ABOVE AND SUMMARY OF FLIGHT RESPONSE 

3 - PRINT ABOVE AND STRESSES 

4 • PRINT ABOVE AND REDUCED VECTOR AND MATRIX 

OUTPUT 

5 • PRINT ABOVE AND FULL STIFFNESS AND MASS 

MATRICES 

ICONTR - ANALYSIS OR DESIGN ROUTINE OPTION (3) 

CHOOSES ROUTINE WITH WHICH VEHICLE WILL BE 
ANALYZED OR DESIGNED UNDER 1CALC*2. 


1 

M 

AW YS 

- SINGLE ANALYSIS PASS 



3 

m 

MWi’ 

- ISOTROPIC 

STRENGTH 

DESIGN 


4 

m 

MWTSTF 

- ISOTROPIC 

STRENGTH 

c 

FLUTTER 

DESIGN 

7 

m 

MWTC 

- COMPOSITE 

STRENGTH 

DESIGN 


8 

m 

MWTCST 

- COMPOSITE 

STRENGTH 

L 

FLUTTER 

DESIGN 

11 

m 

W INPUT 

READ NEW 

INPUT DATA 





(a) Page 1. 


Figure 22,- Inputs for WADES program. 
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UNAL 


IPIOT 


NX 

NY 

^4RAT 
NRLYS 
IFOPT 
I OPT 
ipr 

JfT 

IPO 


JPO 


ANALYSIS CONTROL OPTION VS PLIGHT CCNOtTION ll6*0«A*U 

USI-R MUST INPUT ANALYSIS OPTIONS RiOUlREO FOR 
EACH OF •NOFC* FLIGHT CONDITIONS INOFC<AI I LAST 
FIVE OPTIONS CONTROL SEPARATE ANALYSIS FUNCTIONS. 

FOR J-ltNOFC EXECUTE FOLLOhING CONFUTATION 

ll.O)-l* STIFF - CONFUTE STIFFNESS NATRIX 

LOADS - CONFUTE WORK EQUIVALENT LOADS 
GROSS LIFT*<MaOOY^MTWING«MfUELI*ALOAO 
•2* LOADS • CONSTANT M/S LOADS 

NA^SNX - CONSISTENT NASS NATRIX 
(A.JI*!. NATFR2 - NATURAL FREQUENCY ANALYSIS 

FLUTER - CONFUTE FLUTTER FREQ & MACH NO. 

2f OIVERG - CONFUTE DIVERGENCE Q. 

FCR lANALdtAlt 

U*4i«0f USE ISOTROPIC NATERIAL PROPERTIES 
-1. USE MULTILAYER COMPOSITE PROPERTIES 
•2* USE MULT I THICKNESS FUNCTION CONPQSITE 
I2»4)>|* USE FINITE DIFFERENCE STIFFNESS MATRIX 
(NOT available IN ACSYNT VERSION) 

0.4)«tt OAERO - RECOMPUTE AERODYNAMIC MATRICES 
I4«4)-1. SrOCi - RECOMPUTE T« 0. & C. 

(5f4)*l* WINIT - RECOMPUTE PLANFORM GEOMETRY 
PLOT CONTROL OPTICN (8*0) 

m«Of 00 NOT PLOT APPROPRIATE FIGURE 
(l)*lt PLOT DEPTH FUNCTION DISTRIBUTION 
(2l«l. PLOT THICKNESS FUNCTIONIS) DISTRIBUTION 

(3) >l. PLOT H-OISPLACEMENT VS FLT.CON. 

(4) «1, PLOT FREQUENCY MOOE SHAPES 1-S. 

(5) «l» PLOT FREQUENCY MODE SHAPES 4-6. 

(64-1. PLOT FREQUENCY MODE SHAPES 7-9. 

ONLY UP TO NEIG FREQUENCIES ARE PLOTTED 

(7) -l. PLOT FLUTTER DETERMINANT VS FREQ & MACH NO. 

(8) -l. PLOT PISTON theory SPANMISE SHEAR AND 


MOMENT DIAGRAMS. 

-2. PLOT CONSTANT M/S DIAGRAMS. 

•3f PLOT STRIP THEORY DIAGRAMS. 
number of CHOROMISE DISPLACEMENT FUNCTIONS (3) 

LIMlTSl I TO 6 FOR 'PMLT* 

NUMBER CF SPANMISE DISPLACEMENT FUNCTIONS (5) 

LIMITS* 3 TO 8 FI3R ‘MMLT* 

NUMBER OF STRESS RATIO ITERATIONS (6) 

NUMBER CF COMPOSITE PLYS (0) 

FLUTTER PLOT OPTION FOR X,Y SCALING (2) 

HING PLOT OPTION FOR X.Y SCALING (3) 


NUMBER OF SPANMISE TEAMS IN THICKNESS FUNCTION* (5*3) 
FTT* VS MULTIPLE THICKNESS FUNCTION 
NUMBER OF CHOROMISE TERMS IN THICKNESS FUNCTION* (5*3) 
FTT* VS MULTIPLE THICKNESS FUNCTION 
NUMBER OF SPANMISE TERMS IN FO L FC. (5) 

IF TCR>0« IFD-NUMBER TERMS TO TO USED IN C'MVE 
FIT OF UPPER SURFACE OF WING. 

NUMBER OF CHOROMISE TERMS IN FC C FC. (1) 

IF TCR>0* IFO IS INTERNALLY DEFINED. 

TUTAL number of COFFICIENTS OF FD C FC MUST 


(b) Page 2. 


c igure 22.- Continued. 
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NOVTH 


BE LESS THAN 10. 

> NUMBER Of PLY ANCLE DESIGN VARIABLES lOI 

(LSEO ONLY WITH COMPOSITE ANALYSIS) 

NXSIG > NUMBER Of CHOROHISE STRESS LOCATIONS CHECKED U) 

NYSIG ■ NUMBER Of SPANWISE STRESS LOCATIONS CHECKED U) 

NXGC > NUMBER Of CHORDWISE MINIMUM GAGE CONSTRAINT (S*9I 

LUCATIONS Of THICKNESS CHECKED DURING DESIGN 
NYCC ■ NUMBER Of SPANWISE MINIMUM CAGE CONSTRAINT (5*9) 

LOCATIONS Of THICKNESS CHECKED DURING DESIGN 
NOfC > NUMBER Of PLIGHT CONDITIONS ANALYZED U) 

NwINCS « NUMBER Of WINGS (2) 

THIS EffECTS ONLY THE WEIGHT CALCULATIONS (IE. 

WING WEIGHTS ARE NWINGS*WT)» ALL INPUT WEIGHTS 
(WBOOYtWfUELi ARE BASED ON THE PARTIAL VEHICLE 
WEIGHT involved. 

NCM • NUMBER Of CONCENTRATED MASSES ON WING (0) 

NCMfLT « NUMBER Of CONCENTRATED MASSES THAT MAY BE VARIED (0) 

DURING PLIGHT CONDITIONS 
(NCM*NCMPLT<U ) . (NCHfLT<5) 

NEIC > NUMBER Of EIGENVALUES TO BE COMPUTED (3) 

ITOPT > COMPOS IT'£ «LV TRANSfORMAT ION OPTION (0) 

U • general TRANSf ORMATION 
1 > VECTOR TRANSfORMATION 

ITYPES • MATERIAL TYPE OPTION Of WINC SKIN (0) 

0 - INPUT MATERIAL PROPERTIES 

1 > USE ALUMINUM DEfAULT PROPERTIES 

2 « USE TITANIUM DEfAULT PROPERTIES 

3 > USE GRAPHITE/EPOXY DEfAULT PROPERTIES 

(SEE TABLE 1.1) 

IBUff « ANALYSIS/Or-SIGN VAR (ABLE LINKING fOR MWT^MWTSTf (30... .491 
CEfINES VARIABLE LOCATION IN COMMON/ BWSAV/ TO 
6E TRANSfEREO TO X-VECTOR IN /CNMN2J/ DURING 
DESIGN. OEfAULTS TO THICKNESS VAR I ABLES-f TT 

LEND 


(c) Page 3. 

Figure 22.- Continued. 
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1 BL0CK2 I NAMEUSr CN0ES6N 


I 


CNOESCN 

FTT 

TKiN 

Pti 

TM 

fO 


PC 


TPR 


XTCR 


XTCT 


THEY 


■ INITIAL THICKNESS FUNCTION (L) I1S«0.I 

• MINIMUM GAGE THICKNESS (LI (0.1 

*•* IF (ANAL! U4I>2» USE FOLLOWING REPRESENTATION •** 

• INITIAL MULTILAYEKEO THICKNESS FUNCTIONS (LI S( 10*0. I 

- MINIMUM GAGE THICKNESS OF LAMINATES (LI (0.1 

**• * 


• DEPTH FUNCTION OF DISTANCE OF UPPER SURFACE 

FROM MIOPLANE. IF TCR>0» FO IS COMPUTED 
FROM T/C(RQ0T1 ANO T/CITIPI. (LI 

> CAMBER FUNCTION OF DISTANCE OF CAMBERED 

SURFACE! USED ONLY IN CALCULATION OF PISTON 
THEORY LOADING. (LI 

> THICKNESS FRACTION OF EACH LAMINA WITH RESPECT 
TO TOTAL THICKNESS (TL(II/TCI (USED ONLY FOP 
IANAL(lt4l«ll. 

• X/C LOCATION OF T/C(ROOTI 

IF XTCR>0.; USES XTCR-O.SO (BICONVEX WINGI 
IF XTCR>0.; USES INPUT VALUE IN CURVE FIT OF 
SURFACE 

> X/C LOCATION OF T/C(T1P| 

IF XTCT«0.! USES XTCR FOR ENTIRE WING 
IF XTCT>0. ; USES LINEAR TAPER FROM XTCR TO XTCT 

• COMPOSITE PLY ANGLE ORIENTATIONS IF NOVTH«Ot 

ONLY THE FIRST NPLYS LOCATIONS ARE USED. 


( 10 * 0.1 
1 10*0. I 
( l./NPLYSI 
10 .) 

( 0 .) 

(15*01 


R 

SPAN 

THETl 

THET2 

TCP 

TCT 


LEND 


*** THE FOLLOWING VARIABLES MUST BE READ **** 
** ONLY WHEN USING THE STAND ALONE VERSION ** 


ROOT CHORD (LI 10. I 
SEMI>$PAN OF WING (LI (0.1 
SHEEP OF LEADING EDGE OF WING (OEGI (0. I 
SHEEP OF TRAILING EDGE OF WING (OEGI (0.1 
T/C OF WING AT ROOT (0.) 
T/C OF WING AT TIP 10. I 


IF TCT«Ot ASSUMES CONSTANT T/C VS SPAN 
IF rCT>Of ASSUMES LINEAR TAPER FROM ROOT TO TIP 


(d) Page 4. 

Figure 22.- Continued. 
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I auiCK3 1 namelist 6MATERL 


LMATERL 

RHOS 

RHOC 

MQf 

RVA 

CUTS 

CWTC 

CRAVTV 

XL£ 

XT6 

6XS 

evs 

GXYS 

cxzs 

GVZS 
XVNUS 
YXNUS 
EXC 
6YC 
GXYC 
GXZC 
GY ZC 
XYNUC 
YXNUC 
LENO 


OENSirv OF SKIN MATERIAL < /L*«9l lO.I 

OCNSITV OF CURE MATERIAL t /L*«9I 10.) 

UENSITV OF FUEL (i/L«*3) (0. ) 

FRACTION OF MING VOLUME AVAILAOLE TC CARRY 10.) 

FUEL IN. 

fraction cf non-optimum meicht for skin 10.) 

IMTS • MTS*a«CMTS) 

fraction of non-optimum MEIGHT for core CO. I 

INTC • MTC»I l^CNTC) ) 

GRAVITATIONAL CONSTANT IN USER UNITS IL/SEC**ZI 10.) 
LOCATION OF LEADING EDGE OF STRUCTURE IX/C) CO.) 

LOCATION OF TRAILING EDGE OF STRUCTURE CX/C) CO.) 

CHOROHISE MODULUS OF SKIN lA/L^^Z) CO.) 

SPANMISE MODULUS OF SKIN " 10.) 

INPLANE SHEAR MODULUS OF SKIN * CO.) 

TRANSVERSE SHEAR MODULUS OF SKIN " 10.) 

TRANSVERSE SHEAR MODULUS OF SKIN " CO.) 

POlSSON^S RATIO OF SKIN CO.) 

PCISSON'S RATIO OF SKIN CO.) 

CHOROMISE MODULUS OF CORE CA/L**2) CO.) 

SPANMISE MODULUS OF CORE " CO.) 

INPLANE SHEAR MODULUS OF CORE * CO.) 

TRANSVERSE SHEAR MODULUS OF CORE •• CO.) 

TRANSVERSE SHEAR MODULUS OF CORE " CO.) 

PCISSON'S RATIO OF CORE 10.) 

POISSON'S RATIO OF CORE CO.) 


(e) Page 5. 

Figure 22.- Continued. 
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I 


I 6LQCK4 


STRUCrUAAL DESIGN RLICHT CONDITIONS 


READ THE FOLLOWING STRUCTURAL DESIGN FLIGHT CONDITIONS 
l■l•NOFC FORMAT (9F10.S) 


SSm XMOin Foil I CAMMAIII ALTII) 
ANzm weocviii hfugui) wsm 


WHERE! 

SS - SPEED OF SOUND AT FLIGHT CONDITION IL/SEC) 

XNO > FREE STREAM MACH NUMBER 

PO ■ STATIC PRESSURE IP/L**?) 

GAMMA - RATIO OF SPECIFIC HEATS OF AIR 
alt « HEIGHT ABOVE SEA LEVEL ID 
ANZ ■ ULTIMATE STRUCTURAL LOAD FACTOR 

WBOOV - TOTAL BODY COMPONENT WEIGHTS iLtfS) TO BE CARRIED BY NWiNCSt 
EXCLUDING WING LOCATED WEIGHTSt AND FUEL TO BE LOCATED 
IN WINGS. 

WfUEL • TOTAL WEIGHT OF FUEL TO BE LOCATED IN WING lAI 
wS - W/S* WING LOADING I i/L**2) 

USED AS INPUT only FOR ANAL(2tl)>2 


I BLOCKS 1 CONCENTRATED MASS INFORMATION 


READ THE FOLLOWING CONCENTRATED MASS IHFORMATION 

l>L*NCM FORMAT I3F10.S) 

XCMin YCMIll WCMIl) 


WHERE: 

XC.M • CHOROWISE LOCATION OF CONCENTRATED MASS (X/C) 

YCM > SPANWISE LOCATION OF CONCENTRATED MASS iV/IB/Z)) 
WCM • WEIGHT OF MASS (A) TO 3E LOCATED ON A SINC.LE WING 


I 


READ THE FOLLOWING CONCENTRATED MASS INFORMATION THAT VARIES 
WITH FLIGHT CONDITION 

I-l.NCMFLT FORMAT (5F10.5) 

XCMin YCMIl) AMFLTI I.JIf J-l.NOFC 
WHERE: 

XCM • CHOROWISE LOCATION OF CONCENTRATED MASS (X/C) 

YCM > SPANWISE LOCATION OF CONCENTRATED MASS (Y/(B/2)) 
AfvFLT • WEIGHT OF MASS (i) VS FLIGHT CONDITION TO BE 
LOCATED ON A SINGLE WING. 


(f) Page 6. 

Figure 22.- Continued. 
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I BL0CK6 


NAMELIST CCNSTR 


I 


CCNSTR 

SMART • MAXIMUM TENSILE STRESS IISOTROPIC) li/L**2) (0.) 

WMAX • MAXIMUM ALLOWABLE DISPLACEMENT CONSTRAINT (L) <0.1 

EIGMIN - MINIMUM ALLOWABLE NATURAL FREQUENCY ICPSI 10.) 

FFNIN • MINIMUM ALLOWABLE FLUTTER FREQUENCY <CPS) <0.1 

FMMIN • MINIMUM ALLOWABIE FLUTTER MACH NO. (0.) 

QFMIN • MINIMUM ALLOWABLE FLUTTER DYNAMIC PRESS (i/L**2) 10. 1 
XMAX > MAXIMUM CHOROWlSE SCALE VALUE OF WING PLOT (L) 

XMIN • MINIMUM CHOROW'SE SCALE VALUE OF WING PLOT <LI (0.) 

VMAX • MAXIMUM SPANWISE SCALE VALUE OF WING PLOT <U 

YMIN • MINIMUM SPANWISE SCALE VALUE OF WING PLOT <L) <0. ) 

FOR TRUE SCALE USE ( XMAX«XMINI/( VMAX-YMIN) ■ l.S 
XMAXFM • MAXIMUM MACH NO. OF FLUTTER DETERMINATE PLOT (6.1 

XMINFM • MINIMUM MACH NO. OF FLUTTER DETERMINATE PLOT (1.) 

YMAXFF • MAXIMUM FREQ. OF FLUTTER DETERMINATE PLOT (CPSI (6.1 

YMINFF • MINIMUM FREQ. OF FLUTTER OETERM»NA*« PLOT (CPSI (1.) 

EPSF ■ VARIOUS CONVERGENCE TOLEliANrS» 

1-3 • FLUTTER CONVERGENCE lOLERENCES (3«0.001) 

6 - SCALE FACTOR FOR FLUTTER DETERMINANT (0.5) 

5 > HEIGHT .ONVERCENCE OF MH^ AND MHTC (0.0001) 

6-8 • NOT CURRENTLY U5E0 (3*0. I 


FF • ON INPUTS INITIAL GUESS TO GUTTER FREQUENCY (CPS) 

FM > ON INPUTS INITIAL GUESS TO FLUTTER MACH NUMBER 

(END 


I BLOCKT I NAMELIST &CNMN 


(READ BY C0MP30) I 


CCNMN 


IPAINT 

NOV 

m 

PRINT CONTROL OPTION (NO«0»YES>0) 
NUMBER OF DESIGN VARIABLES 

(0) 

UMAX 

• 

MAXIMUM NUMBER OF ITERATIONS 

( 10) 

NSlOE 

m 

SIDE CONSTRAINT OPTION 

(0) 

NSCAl 

m 

SCALING CONTROL PARAMETER 

(0) 

NFOG 

« 

GRADIENT CALCULATION CONTROL 

(0) 

FOCM 

m 

RELATIVE FINITE DIFFERENCE STEP 

(0.01) 

FOCMM 

m 

MINIMUM RELATIVE FINITE DIFFERENCE STEP 

(0.01) 

CT 

m 

CONSTRAINT THICKNESS PARAMETER 

(0. 10) 

CTMIN 

m 

MINIMUM ABSOLUTE VALUE OF CT 

(0.006) 

CTL 

m 

LINEAR-CONSTRAINT THICKNESS PARAMETER 

(0.01) 

CTLMIN 

M 

MINIMUM ABSOLUTE VALUE OF CTL 

(0.001) 

TM6TA 

m 

PUSH-OFF FACTOR 

(1.00) 

PHI 

9 

PARTICIPATION COEFFICIENT 

(5.00) 

OELFUN 

9 

RELATIVE CONVERGENCE OF OBJECTIVE 

(0.001) 

OABfUN 

m 

ABSOLUTE CONVERGENCE OF OBJECTIVE 

(0.001) 

URM 

m 

CONSECUTIVE ITERATIONS OF REQUIRED CONVERGENCE 

(3) 

VLB 

9 

VECTOR OF LOWER BOUNDS ON 'X* 

(30*0.) 

VUB 

9 

VECTOR CF UPPER BOUNDS ON *X* 

VUB C VLB ARE USED ONLY FOR NS10E>0 

(30*0.) 

SCAL 

9 

VECTOR OF SCALING VALUES FOR «X* (X-X/SCAL) 
USED ONLY FOR NSCAL<0 

(30*0.) 


LEND 


(g) Page 7. 


Figure 22.- Continued. 
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I BLOCKS I COMPOS ITi PLY NATEPUL PPOPEPTIES (BEAD BY COMP2S) I 


PEAO THE POLLOWING COMPOSITE MATERIAL PROPERTIES POR 
1ANALIU4I > 0 AMO lTVPES*Ot PORMAT ISP10.3I 

RE AO MODULI AMO POISSON'S RATIOS 

ELll EL22 NUL12 NUL21 GL12 

GL13 CL23 

READ LAMINA DENSITY* PLY THICKNESS. VOLUME RATIO* PIBER CONTENT 

RHOL TPLY KV KP 

"EAO LIMIT STRESS/STRAIN CONSTANTS 

SLUT SLllC SL22T SL22C SL12S 

SL23S PHIOEL 

MHERE THE VARIABLE ENDING IN *T> OR 'C DENOTES TENSION OR 
COMPRESSION IN THE APPROPRIATE COMPONENT DIRECTIONS AND *S* INDICATES 
SHEAR. POR PURPOSES OP DESIGN SL13S-SL23S. PHIDEL IS THE STRAIN 
ALLOWABLE POR DEL AMINAT ION* WHICH IS NOT USED I*‘ THIS PROGRAM VERSION. 


I BLKK9 I PLY ORIENTATION DATA (READ BY COMP14) I 


IP NOVTH>0* SKIP INPUT PLY ORIENTATIONS IN NAMELIST CNOESGN 
IP N0VTH>0* READ COMPOSITE ANALVSI3/0ESICN VARIABLE 
TRANSFORMATION ACCORDING TO ITOPT. TRANSFORMATION LINKS 
PLY ORIENTATION AND THE DESIGN VARIABLES ACCORDING TO{ 
ITOPT-0, THETiiJ • SUM T(I*J)*XU) ♦ THETII^NPLYS) 
ITCPT-l* ThET(I) • T(I)*X4I) ♦ THETfUNPLYS) 


FOR IT0PT»0* 

HEAD J-ltNOVTH 
T< I*Jt *1*1*NPLYS 

POR ITOPT-l* 

READ T( II* I-l.NPLYS 

r.EAO CONSTANT PORTION OF ANALYSIS VARIABLES 
THET(I*NPLVS) *1*1*NPLVS 
READ INITIAL DESIGN POR THET 
THETII ).I-l*NOVTH 


(h) Page 8. 

Figure 22.- Continued. 
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SAMPLE DATA POA P-SA INPUT TO ACSVNT 


P-SA/B ACSYNT-MAOES NINC HI ICHT CORRELATION 
AOPTNS IPRT>l*IANAL>ltlfOtOtOf lti*0»0«0t 5*0* OtA«U 
IPLOT*li Itlt0i0.0.0*2i IPRNT>l.JOUMP«3tlCONTA«3t 
1FT>A*<»*0*JPT«A,4»0« IF0-S*NCN«0*NCNPLT>2»N0PC>2t ITVPES*lf tlNO 
INOESCN FTT«0.4f0.fU.t0.*0.Atl0*0.« TNIN«0.02i UNO 
AMATIRL RHOP«O.0S6»PVA«O.0»CRAVT .'■JSo.JOTf 
XLE«O.19tXTE«0.S5iCHT$>0.0»CHTC>0.0»XVNUC>0.0iVXNUC«O.0* 
EXC>A29000.»E VC* 77BO0O.»GXVC*0.tCX2C«394O0.«GYZC* 116600.. 
RHOC*0. 00666* 6EN0 


12360.0 

1.09 

6.33 

1.6 

21900. 

9.6 

9997.0 

0.0 

0.0 


13392.0 

0.22 

16.? 

1.6 

0.0 

6.0 

10632.0 

0.0 

0.0 


0.59 

0.60 

266.0 

-6327.0 

LAN0IN6 GEAR 

O.IO 

1.0 

206.0 

119.0 

AIN-9B/TIP TANK 


•CNSTR TCACE«0.02.WMAX*30..E1GM|N-3.0*PFNIN*3.0.PMMIN>2.50.0PMIN>23. 2. 
FN>3*3.0tFF> 3*12.0. XMAXFMa6.0«XNINFM*1.0.YMAXFF>16.U.YMlNFP>6.0. SEND 
ftCNMN IPP|NT-0.NSCAL«0. ITMAX*30.NFOG*O.NSIOE-0. »EN0 


(i) Page 9. 

Figure 22.- Concluded. 
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original PAQt JS 

OF POOR QUALITY 


F-5A/A »iiNG F.EIGHT cn«Het«rtuN • ne^AULr stamtimg ftt 

iOl»T>‘iS IP«T»l» IP«NT*|» JDiJNp.4, ICLtHTAaJt 
I fi| 1 « ) *o « lit II tC * 2 1 

MPLVb»'i*IP0«*i*'.C**«O*<MPL1«A«f.UPC»3f ITyPeS-1 1 ‘f'*0 
♦ mOCSGm «•! 3*.itSP#^ •I»I •9tT.«6n«^2.0»T«eT2««5*0f7CP«0«f»<**t 
FTT«»J*e*»j»0«^’*0*<'*‘>*U*8* lM*0»0f T'<|N«0»02f ♦€<'•0 

i**AtfePt «HuC««i«OOOiPiy«MHtif ■/l•056•PvA■O•0tGRAvrY■^e6•307«XLe■0•15*X7e■0•9S• 
t«C»*29OH(>,tfcrc»77Hn0O**rf«VC«rt,O,G*IC*J^^'7‘**»GVlC*n»ft')0, ♦ « HOC •0, 00848 1 
xym'C«o»o»Y*MOC«o*ot *6^’o 


12360,0 

l.'»3 

6,33 

1.6 

21500 , 

0.8 

95 ^ 7 , 

« l,0 

0 ,(* 


13392.0 

0.22 

14,7 

1.6 

0.0 

<►.0 

l <><* 32.0 

0.0 

1.0 


13392 . i > 


\ 6,7 

1.6 

0,0 

9.8 

9359 , 

0.0 

0,0 


0.55 

0,40 

748 , 

• 4327 . 

24H . 

O.lo 

1.0 

204,0 

115,0 

U5.0 

0.3 

0,56 

0,0 

0,0 

112.75 

0.3 

0,75 

0,0 


119. 0 


»Ca. 5TR T&AC.^*0.0?fi«.'»&»«lOO,F 

HGMlN - j . OfPf - , PMmin « 2 , 90 * i.'FM '«» 23 . 2 » PM « 3 * 3 . nO » PP « 3 * l2 .'). 

A -llNpHol ,0f yM&AP»-«18*0,yMl(<pp»l».0* 

iPHt^rNii, r r**Ax«30t »tft' 


Figure 23. • Seunple of input card to WADES program. 
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li « 0 I • r***cxwio*c •f*'* i«^wr • «c«»iv«t i/t/r» 

alMft aClt*** eO*«|k«TtP*> • P|*«Ulf ff1 


lh0tiT COir»^^ •*|n»«CTI*t 


I 


ttiA’.TlII COMTtOk O^TlOtt lt«*»*kWMe««,VClM) 
•Tu» kO«a* 

I I I I 0 

I t u u • 

I t 0 V 4 

eu**i*<*i ••• B 4 fc»o iToc MO" 

0 I t t I 


*kOT CO"T«U 0 »TfwM (lOcOrMNO« 0 **CI»«) 

9 <*T*i tmicr r^kC’" 'MOi rM«* 

0 I t U 4 

'Mmi nr OCT kO«oi 
0 0 • 


0k«C4I*f|nM (0C4*UkT) (VaR) V 

4kUl 




Of CMORO"t»< •Ok»"k"l4tl 

<»» 

(NS) 





"0"UC» 0» "OkTMO"i«k| 

«k> 

(nV) 





.ir STMM 44TT0 |TC*AriP«.» 

<•) 

(NRaTT 





MO"»4* Of CO""00|U "wtl 


(N^kvC) 





IfV SCALING 0P1 

»>» 

(lotn 





*.» •C4tl;<« ORTlON-»tuTT»K OCT, 

(<) 

ri4tf4T) 





NU"Ui« Of fit r|»M» • 4T|.(1) 


(14T) 

i • 

s 

0 

0 

'•U-**iC<4 Of II T|4m4 • tTkltl 

(!) 

(J4T) 

i 0 

0 

0 

p 

ToT4t Of Tt»"l • Mk(I) 


<"»t ) 

• 0 

0 

p 

p 

•,U"ut<< or |T4 T|»n» > fO,fC 


(14U) 





*tU"UC<« 04 I| T{u,.| • >0,4C 


<JSO) 





TqTiL **0»*mR 0* Tt*"l • 40.4C 


(««rO) 





0C6"kkK 04 4MK4Dr’< 04 ' 


(NK) 





•*0"n*" of C"0MO«l4i » *•,•» kOC* 

(i) 

(NVlTG) 





tOi^bk- 04 |44««ft{|t kPC. 

(Ml 

(N>SIQ) 





Nu"viii jf e"0*OotaC OkO'^tTaic lo>* 


(••*nc) 

% 0 

0 

0 

p 

**U"»C^ 34 4"4h«li|A i*|o*CtMlc k9c« 

(S) 

(••^GC) 

% 0 

0 

0 

0 

"U""!" 'of -IN(;« 

(4) (MMtnnC) 





•«0"»t“ <.f »kU"T C0»<0IT10**$ 

(n 

(N04C) 





04 C0*<Ce*«T-4t4C ma«S(1 


(NC«) 





•*y"b»" j4 cO**C. "4»*k* / 4klC».T fO", 

(»(C"RkT5 

4* 




<«li"kC" of 41<;c*«v«cf U«« "Ctiiaivp 

(!) 

C*«tt6) 





'4u'*tlk>« 04 4k* 4<»r.kk OtSIC'* V44«A4kCt 

COvtm) 





ANAL^d! VAMlAdtt Ll^P!*<0 


riC044) ] 

iu II 

IS 

11 

IS 


(mi rmgm 1. 

rigurm 24. ~ Input oopy gunnrntnd by ooaputnr. 
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WIOINAL PAGE IS 
^ POOR QUALITY 


• *T>t« t^*ul • l/|/t» • 


*k««**0** v«»t4%wtK HtOCX 3 

*viuf c«*0«0 ■ tl«.% • IftltV 

ki4tft*i4 kuei *»lk» ■ U.tf iour 




rNtcy«r»0 « icofM 

Cofo* iClIHTk 


0TT 

■ 





•000l-'<.' u. 

A. 

0. 



OOAOk aO A. 



Tm 

■ 

lOAOfaOl 



n 

• 


0. 

0. 


I* 

C. 

A. 

Ot 



0. 



fC 

> (»• 


0. 

0. 



0. 

0. 

0. 


0. 

0. 





4W^k**»4H •ca«444i«T|f{ON 



^/C(4nuT) • 

T/W(ri»*) ■ 

w.OOOOW# 



• O.PftAOOO 

4*oOC r/C(T) a 

V.UOOOOO 


uCNiITtCI (rcl0toT/L(*t6T»>«*}) 


CC^t *»!•< 

.oo««#o .tooooo 

VUwUMf 4W*Tk4tC» • 

«lUN*n*T I •l»« nflimt, |k{s - 
««4V|r» H.»’<0TN/|CC»*i) 
10C4T|W«* n» ,rNR ROy. t.k,* 


»w»t 

, 0 $* 0 Ab 

4, AOOU 

1«4R0f CO**! ■ O.RA^k 

• l*«0 »,f ,• (t^UO 


■LOCK 9 


Ik 

It 

01 * 

6 it 

0 v 2 


** 0 C Jtw* 

»• 


• 4 T| 0 Ut 

co»e 

«.«^00C*t/9 

».T0ttOl.OS 

, 4 **nk**<> 

.na 0 k*w» 


t .o«ooi*Ar 
. 4004l*4T 
.»<tonr »f7 

• *0wi>l*4 t 


•*0U90N's 4«riO0 
CO*t 0«iM 

NU«y < 1.000000 (Soooeo 

NIITI O.OOOOOO .100000 


»klO“T CONRlTIONl . N04f a 


I^CiO 0^ *C"A,o 


MkCw s.i(««4C<» 


1 .nuOu 

4t. »l*i.'Ot 

• Oo^ -tlv#"! 

• f • w 

•k0wl»k*» fito 

S^ttO 0^ 4no.«i> 

1 

»4ra a'!a«k" 

4 , 100 4 

1 • 4 0 V g 

4tTii •<:>'' 

• OO^ •€ 


4cOuIat4 yi'tii. 

SntO SOVNO 

m«^tO 


\«A **^4 

t «g000 

4tT|Tl.0l 

•Oo^ •t!c-r 


•kSik<l0lo »"fl. 


^ ILOCX 4 


1 • U^gO 

static 

OtilOO 

iMfiU. 

LJaO ^4CI0« 

OtOOOO 

^.0 
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APPENDIX A 


CALCULATION OF F-5 EQUIVALENT CORE PROPERTIES 

Th« rib and sp^r matcriml and structural propsrtlss sho%m in ths fol- 
lowing tables wars used to csiculats ths squivalsnt core properties in the 
F-5 correlation. The values selected here are only representative of actual 
section properties, since the actual substructure varies with location on 
the wing. The general approach taken in selecting these properties was 
either to pick a location which was typical of the rtiffness or mass in a 
critical region of the wing or to obtain an average value of the particular 
property if applicable. For example, the average density of the ribs and 
its equivalent distributed value could be calculated directly, since the 
weighs of each rib was known. The general assumptions used in deriving 
the equivalent distributed section properties are outlined in the section 
on equivalent core propert-iea. Taioles A-I and A-Il contain a summary of 
the representative F-5 spar and rib section properties used in these 
calculations . 

The equivalent core shear stiffness in the y-direction was computed 
as a component of the shear stiffness along the structural reference axis, 
x/c “ 0. 35. The actual shear stiffness of the core is approximated as the 
product of the shear modulus with the sum of the material areas of the F-5 
spar shear webs at WS 64 : 

“y. - 0E\,,b 

- (4.0X10*psi) (0. 562 0.208 f 0.288 

f 0.428 ^ 0.759 ^ 0.324)in* 

- 11. 76X10* lbs (A-1) 

The equivalent distributed cross-sectional area at WS 64 is: 


110 


\c " 


e 

S 


- (25.78 in) (3.6 in)cos 20.5° 


- 06.9 in* 


(A- 2) 


Aft«r rotation to th« structural rafaranca lina, tha aqu'^valant cora shaar 
modulus in tha y-diraction bacomas: 


0 


yxc 



CO8*0^ 


_ (11. 76X10* lbs) (COB 20.5°)* 

06.9 in* 

- 110,600 psi (A-3) 

Tha aquivalant cora shaar stiffnass in tha x-diraction was computad 
as tha sum of tha shaar stiffnassas of tha individual ribs at (x/c) * 0.4. 
All ribs wars aluminum axcapt tha landing gaar rib &t KS 64, which was 
stsal. Tha shaar stiffnass of tha P-5 ribs is: 

<°^x)p5 - ®\.ab 


- (4.0X10°psi) (1. 329 0.63 -t- 0.216 ■*- 0.22 ^ 0.84)in* 

- 1.935X10'^ lbs (A-4) 

The equivalent cross-sectional area of structural core in the x-direction 
is : 


^c • <^ip ^ ^root)^*/^ 


- (1.42 in 5. 5 in) (142 in) 


- 49i. 3 in* 


(A- 5) 
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Th« •qulvAlant cora shaar modulus in tha x-directlon bacomas: 


(OA„), 


- (1.935Xlo‘^lba)/491. 3 in* 


- 39,400 pal 


(A-6) 


Tha aqulvslant cora banding stiffnasa in tha y>diraction waa computad 
as a componant of tha banding stiffnase along tha structural rafaranca 
axis. Tha structural banding stiffnasa of tha cora is computad as tha 
product of modulus of alasticity with tha sum of tha momants of inartia 
of tha P-5 spars: 

<="x)p5 - = (^x>spars 

- ( 10. 4X10® psi) (1.884 0.867 + 0.867 

+ 0.867 + 1.343 + 2.233 + 1.062) in*' 


- 9. 488X10'^ lbs in* 


(A- 7) 


The equivalent distributed moment of inartia of tha core at WS 64 is : 

^xc “ 0-3c(t^^)*(cos 0,)/12 

- (25.78 in) (3.6 in)’(cos 20.5°)/12 


- 93.88 in* 


(A-8) 


After rotation to tha structural rafaranca line, tha equivalent banding 
modulus in the y-diraction bacomas: 


*“x>F5 . 

co.‘ 9, 

xc 


- (9. 488X10“^ lbs in*) (cos 20. 5)*/93.88 in* 


- 777,900 psi 


(A- 9) 


aqulvalMt cor* bonding *tlffn*** in th* x-d.ir*ction w** conputod 
*■ th* auBi of th* bonding otlffn*** of th* rib* *nd * componont of th* 
bonding stlffn*** of th* apar* along th* otructural axl*. All rib* war* 
alunlnun oxcopt th* landing gaar rib which wa* ataal. Th* rib contribution 
to cor* bonding atlffn*** 1* conputod a* th* aun of th* products of th* 
momonta of Inartla of th* P»5 rib* with tholr moduli of oloaatlclty: 




<®Vrib* 


- (10.4XlO®p*i) (5.69 4 5.28 t- 4.79 -► 0.997 4 1.504 
4- 0.509) in* 4 (25Xl0®p*i) (5.28 in*) 


- 2. 722X10* lb* in* 


(A- 10) 


Th* oqulvalont dlstrlbutod cro**-**ctlon moinant of Inartla for a trap*- 
zoldal aroa Intogratod from th* root to th* tip la: 


I - 657.4 in* (A-11) 

yc 

Th* oqulvalont banding modulus of th* cor* in th* x-dlr*ctlon bocomo* th* 
sum of th* componants from th* ribs and th* spars: 


xc 


(El ) 5 (El ) 5 

— JLX2 ^ — j2L-L2 

yc ^xc 


- (2. 722X10*lbs in*)/657.4 in* + 15,200 psi 

- 429,400 psl (A- 12) 

Th* oqulvalont coi* donslty of th* F-5 was computad as th* sum of th* 
donsltlos of th* rib* and th* spars. Th* donslty of th* rlbe was computad 
from th* waight of th* ribs dlstrlbutod ovar th* volum* of th* structural 
planfom. Th* donslty of th* spars was computad from th* croas-sactlonal 
aroa of th* matarial in th* spars at WS 64 dlstrlbutod ovar th* cross- 
sactlonal aroa of th* wine structur*. All spars ar* consldorod to b* mad* 
of aluminum, "0.10 Ibs/ln^. Th* aroa of spar matarial at WS 64 Is: 
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- 0.940 * 0.460 0.460 0.460 ••• 0.678 1.259 -f 0.544 


•par 


- 4.801 in* 


(A- 13) 


Then th« density bscomas: 


•pars al A. 




ye 


- (0.10 lb«/in*) (4.801/86.9) 


- 0.00552 lb*/in* 


(A- 14) 


Tha astiaatad voluma containad in tha structural planform froa the WADES 
prograa is: 


''eor. ■ i”’ 


(A-15) 


The total of F-5 rib weights is: 


"rib. * 


(A- 16) 


Tha equivalent density of ribs is than: 


^ribs " '^ribs'^cora 


-80.6 lba/27,200 in* 


- 0.00296 lbs /in* 


(A- 17) 


Thus, tha net equivalent density of tha core for tha F-5 becomes: 

^c " ^ribs ^spars 


- 0.00848 lbs /in* 


tA-18) 
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APPENDIX B 


CALCULATION OP P-S NABOIM8 OP 8APBTY AMD SUMPACE PIT 
OP UPPER WmO SKZN-THICKME88 DISTRIBUTION 

'Hi* analysis dsscribsd bslov was ussd to calculats ths sr'’facs fit of 
ths uppsr wing skln-thicknsss distribution and to calculats ths margins of 
safety for ths critical loadings spscifisd in rafsrsnes 1. Ths calculations 
and results of this section %rsrs ussd to generate the equal skin- thickness 
contour plots of the surface fits in figures 3 and 4. The points used in 
this analysis contain only a representative portion of the total number of 
panels examined by Northrop. 

Table B-I is a sunsuury of the input used in the following calculations. 
The first three columns contain the panel number and the nondiMnsional 
chordwise and spanwise locations on the wing surface. The above locations 
are ratioed to the root chord and the semispan respectively. In the fourth 
column is the Northrop flight condition number identifying the critical 
load condition used to size the particular panel. The parameters B and T 
are the local panel width and thickness. In all cases the length of the 
panel will be considered along with respect to the width for the calculation 
of the buckling stress. The columns PCL and PS give the critical edgewise 
stresses on each of t^ie panels obtained from examination of all the F-5 
loading conditions. The negative sign on the values of edgewise stress 
represents compression. The last column, TCAP, gives the average total 
skin thickness at the adjacent spar locations. This value includes the 
wing skin thicicness, the spar cap material, and the skin landing material. 

The basic strength design and check of the margins of safety of the 
F-5A/B wing skin were carried out on the basis of a combined buckling fail- 
ure analysis. The wing skin was idealized as a series of individual panels 
bounded by the spars and rib*». Bach panel is assumed to be a long, flat 
rectangular plane of unilorm thickness with simply supported edges loaded 
in shear and compression. Allo%#able panel buckling stresses have been 
determined using center panel thicknesses and panel widths between rivet 
lines. The effects of taper in panel width and thickness were approximated 
L;' checking stresses at each end. 
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ITi* allowAbl** us«d and Mrgins of oafoty woro obtainod froai tho 
following oquationa: 


tttMM - gwgytttign 



n 



( 1 - 1 ) 


whara 


B - 10.«ri0*pai 
K - 3.62 


Buckling atraaa - Shaar 


whara 

4.9 i K, i 5.75 

baaad on b/a ratio. 

Tha allowablaa P^r Pa^^, (tha compraaaiva and ahaar allowabla 
atraaaaa), ara obtainad from intarp>olation in figuraa B-1 and B-2, Pgjp 
varaua Pcr^^ ^*cr varaua Pacr/^* 

Tha margina of aafaty ara computad from tha intaraction aquation for 
a panal aubjactad to combinad adga ahaar ar.l compraaaiva atraaa (raf. 1) : 

■ ^•'5 (B~3) 


whara 


'cl 


i£i- 

^max 


and 
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Th« nAX'gin of sAfoty bocoM'i 



Tabi* B-ZI prooont « suaMry of tho calcul«tion of tho allowablo 
•trossoa and aarglna of safaty for a nuabar of rapraaaiitativa panals on 
tha uppar akin of tha P-S wing. Tha valuaa of tha scrasaai^ and aargina 
of aafaty aaaantially duplicata tha valuaa shorn*, in tha Northrop P>t wing 
analysis (raf. 1). 

To obtain diract coaparison with tha f-S matarial distribution, a 
laast^squaras fit of a poiynoaial function in tha spatial coordinatas 
( and was aada of tha uppar skin thicknass and of tha affactiva skin 
thicknass as dafinad in aquation (B>1). Tha polynoaiial coafficiants fron 
thaaa functional fits wars than usad in a datailaU analysis of tha P>S 
by tha MADB8 program. Tha thicknass function usad to approximata tha 
matarial distribution was of tha form 

tij(^n) - (b-5) 

vrhara tha function w^ is dafinad by 

P 

Wp(?,n) • (t ” tan ♦ tan f T 

For tha P-SA/B wing configuration, tha following valuas wars usad: 

R - R - 134.5", S - 151.5", 9 - 32®, and S - -5® 

Tha following analysis two functions wars usad in ordar to avaluata 

tha sansitivity to tha numbar of fraa coafficiants nacassary to astablis.i 
a good fit. Tha two functions ware : 

PT - C f C L -t- 

33 1 4'» s'* 

♦ c 1 + c 

2 ■ 

* C *1* (B-7) 

3 
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and 


FT. 


44 


• c ♦ c 






10 ^ 


4 c^n* 4 




(B-8) 


Th« sacond lunction was also usad to computa tha surfaca fit of tha affac- 
tiva skin-thicknass distribution. 


Iba laast-squaras functional fit of tha abova aquations Is obtalnad 
by solution of tha walqhtad sat of aquations which mlnlalsa tha rasldual 
arror : 


Zyf^ 

4 


t If C 4 
I *■ I 


t If f 

a a 1 


t If f„ 4 t If f 


1 M 


a M 



(B-9) 


Kara tha si*anatlon Is cosiputad for tha numbar of obsarvatlons and tha f)^ 
ara tha spaclal functions corraspondlng to tha k'th fraa coafflclant of 




(B-10) 


Tha function y Is computad from tha actural F-5 skln-thlcknass valuas: 


y • t 


^min 


(B-ll) 


Valuas of t corraspondlng to tha actual F-5 panal thicknassas wara usad 
for tha initial surfaca fits. For tha surfaca fit of tha affactiva thick- 
nass distribution, valuas of t ■ wars usad. Tha systam of aqua- 

tions (B-9) is solvad as a llnaar systam by LDU dacomposltion. 

Table B-IIX Is a summary of tha results obtained from tha laast-squaras 
surfaca fit of functions described in aquations (B-5) through (B-6). Tha 
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co«fflcl«nts hav« b««n conputsd for a nominal valua of ■ 0.02. Tha 

tabuiatad paramatars T33 and T44 ara tha computad valuaa of t corra- 
■ponding to aquationa (B-7) and (B-8). Tha column PTE contalna tha 
computad valuaa of t^ff ualng aquation (B-8) to ganarata tha surfaca 
fit. Tha column SIO(VM) la tha von Mlsas* atraas aa computad from 
fg]^ and fg obtalnad from rafaranca 1. 
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Of POOR 


TABLE B-Z 

TABLB or CBZTICAL DBBIOV PARANBTBB8 


X/C 6fA FlT con 0 T FCL FS TCAP 


1 

0.290 

U • ^ J 

123C-5 

6.27 

C.208 

-40941. 

-53. 

C.357 

2 

0. 340 

0« u J 

123C-5 

6.49 

W.232 

-48261 . 

-293. 

0.395 

i 

0. 390 

0 • u J 

123C-5 

6.32 

0.245 

-52103. 

-363. 

0.429 

4 

0.440 

c *00 

123C-5 

5.75 

0.258 

-54817. 

-2*>»80. 

C.S08 

4 

0. 490 

u* Uw 

123C-5 

5.51 

0.288 

-566 94. 

-2164. 

0.538 

0 

0. 240 

0 • 1 

12JC-5 

6.53 

0.2 42 

-46184. 

7872. 

C.429 

7 

C . 300 

O.l^ 

123C-5 

6.6 0 

0.255 

-50745. 

7918. 

C .46 7 

6 

C. 360 

w. 19 

123C-5 

5. 77 

0.258 

-55530. 

6930 . 

0 . 50 8 

9 

0.415 

0.19 

123C-5 

5.51 

0.288 

-57037. 

3118. 

C . 56 8 

10 

0. 240 


104 

4.73 

0.185 

-45263. 

“5643 . 

0.291 

11 

0. 3C0 

0.45 

1 J4 

4.86 

0. 205 

-44761 . 

-6054 . 

0.321 

U 

U.360 

0 ^ 

104 

4.17 

169 

-47875. 

“7680. 

0.325 

13 

0. 415 

u • V 5 

1 04 

4.15 

0.182 

“4 82 48 . 

-55o8 . 

0.355 

14 

J.240 

0.‘>9 

104 

3.96 

0. 1 50 

-35706. 

-3660. 

0.257 

15 

0.3o0 

0 .59 

U4 

4.08 

0.165 

-37005. 

-4915. 

0. 173 

lo 

J.3O0 

0.59 

104 

3.31 

0.133 

-43701. 

“06 7 7. 

0.287 

17 

0.415 

0 .59 

104 

3.44 

0. 145 

-43214. 

-6 5 45. 

0.327 

Id 

0 . 550 

0.59 

382 

7.50 

u.l 10 

-1287. 

6iI27. 

0.3C8 

19 

C. 180 

0. 72 

361 

3.19 

0.125 

-262 75. 

-04 96 . 

0.250 

20 

0.240 

U. 72 

361 

3.32 

0.130 

-30720. 

-8325. 

C.250 

21 

C. 300 

0. 72 

361 

3.50 

0.130 

-34304. 

-4269. 

4.250 

22 

0. 360 

0. 72 

3ol 

2.62 

0.135 

-34894. 

-9401. 

0.252 

23 

G.415 

0.72 

3ol 

2.90 

0.135 

-36065. 

-8995. 

0.289 

24 

0. 195 

0. 

3ol 

3.14 

0.128 

-15807. 

13523. 

0.270 

2* 

0. 285 

U. dO 

iol 

3.71 

0.132 

-186 34. 

Idl21 . 

0.270 

2*s 

0.375 

0 .od 

361 

3.94 

0.130 

- 1 96 46. 

16792. 

C.27C 

27 

0.475 

0. o d 

3ul 

3.14 

0.124 

“1 434 7 . 

19 3o5. 

0.270 

28 

0.545 

O.dd 

3ol 

3.51 

0. 113 

-i«.n3. 

19947. 

0.270 

29 

0.620 

L . O 3 

361 

3.68 

0. 1 12 

-13362. 

17426. 

0.295 
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TiABLI B-XZ 

CQHPUrSD CBITICXL 8TKBBBB8 AMD NABOim OT 8APBTY 


1 

FCR/E 

FCH 

FNAX 

FitR/E 

FSCR 

RCL 

R$ 

M.8. 

1 

41830. 

41830. 

49549. 

66443 . 

38631. 

0.826 

0 .00 1 

C.210 

2 

485 72. 

48572. 

51071. 

77151. 

40058. 

0.945 

0.JU7 

0.058 

3 

53670. 

92936. 

53282. 

85250. 

41030. 

0.978 

0.CC9 

0.023 

4 

76525. 

63 957. 

63957. 

121552. 

44724. 

0.857 

0.06/ 

C.160 

5 

103844. 

69000. 

89000. 

164945. 

46000. 

0.822 

0.04 7 

0.213 

h 

32204. 

51763. 

52490. 

82920. 

40750. 

0.880 

C. 193 

0. 086 

1 

56740. 

55114. 

55114. 

90126. 

4lol5. 

0.921 

0.190 

C.043 

0 

75995. 

03794. 

63799. 

12J710. 

446 5 7. 

0.870 

C.155 

C.llS 

9 

103844. 

69000. 

69000 . 

lo4945 . 

46000. 

0.827 

C. C68 

C.2U2 

10 

58146. 

5oJl 3. 

56013. 

92 359. 

41883. 

0 . 8o 8 

0.135 

0.205 

11 

67629. 

o0804. 

60804. 

107422. 

43442. 

0. 736 

0.153 

0.3o4 

12 

62431. 

5o4l 5. 

58415 . 

99165. 

42616. 

0.020 

0.180 

0.166 

13 

73105. 

62 856. 

62856. 

116119. 

44290. 

0. 768 

0.126 

0.269 

U 

5 45 3 7 . 

53629. 

53750. 

86626. 

41195. 

0. 664 

0 . C 94 

0.477 

l5 

62165. 

58282. 

58282 . 

98743. 

42574. 

0.635 

0.115 

C.526 

It 

6 L 3 6 b . 

5 7884. 

57884. 

97477. 

42448. 

0.755 

0.157 

0.272 

17 

6 7533. 

60 764. 

60764 . 

107269. 

4342 7. 

0.711 

0.151 

0.348 

la 

81 76. 

ol 7o . 

2 82 1 7 . 

12987. 

12987. 

0 .046 

0.687 

0.4u 7 

19 

58363. 

5o 152. 

56152. 

92 703. 

41924. 

0.4od 

0.155 

0.943 

2C 

58278. 

5o098. 

56098. 

92569. 

41908. 

0 . 5<.8 

0. 199 

0.634 

21 

524 38. 

91951. 

52617. 

83293. 

40795. 

0 .652 

0.227 

C.382 

22 

1009 17. 

o9000. 

69000. 

160296. 

46000. 

0.506 

0.204 

C.730 

23 

82370. 

o571 1 . 

65711 . 

130836. 

4545o. 

0.549 

0. 198 

C.632 

2<* 

631e2. 

587el. 

58781 . 

100327. 

42733. 

0.269 

0.316 

1.091 

25 

48117. 

48117. 

50999. 

75429. 

39971. 

0. 365 

C.403 

C.599 

26 

41380. 

41380. 

49414 . 

85728. 

38516. 

0.398 

0.4 36 

C.475 

27 

59276. 

56737. 

56737. 

94154. 

42099. 

0.2.->3 

0.460 

J.657 

2b 

39395. 

39395. 

48782. 

62575. 

38012. 

0.293 

3.52 5 

C.446 

29 

35208. 

35208. 

47275 . 

55924. 

36866. 

0.283 

0.473 

0.575 
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ORIGINAL PAGE IS 
OF POOR QUALITY 





TABLE B-ZZI 





siBBiABY or BimrACB rZTB 



IFT - 
l.«281 1 
-4o. 1035 

iU4FACE FII OF 
4 JFT • 4 

I l.4i5d 
•*. 10^87 

rnlCKNESS FUNCTION 
NFT ■ 10 THIN > 0.020 

-0.52732** 33.8763 

27.2395 1.29676 

-3.04C35 

•2 5.4222 



SUHFACt FIT UF 

thickness function 



IFT - 

3 JFT • 3 

NFT > o THIN « 0.920 



1.36769 

0.99/52 

7.630*72 -2.38482 

-15.9313 

3.(1314 


SUKhACt FIT uF 

EFFECTIVE SKIN THICKNESS 1 

FUNCT ION 


IFT . 

4 JF T ■ <» 

NFT ■ 10 thin ■ 0.020 



O.7)u53o 

7.2S953 

•14.4\^40 39.5250 

3. 75598 

4.35044 

-2 7, 4‘#64 

- .42 7o 

-11.1797 18.5278 





TAbCt UF 

THICKNESS FIT 

VALUES 





I 

x/C 

blA 

X 1 

T 

T33 

T44 

5lC( VM) 

TtFF 

FTE 

1 

0.290 

O.JOJ 

0.290 

0.2080 

0. 2218 

0.2096 

<t094l. 

J.2318 

0.2317 

2 

0.3 40 

O.OOJ 

0 . 340 

0.2320 

0.2337 

0.2256 

48262. 

0.2571 

0.2528 

3 

0. 390 

0# 000 

0.39w 

0.245U 

0.2548 

0.2431 

52104. 

U.2732 

0.2745 

4 

0 . 4*»0 

U« JJJ 

0.440 

0.2580 

0.2708 

0.2635 

54698. 

0.3015 

0.3004 

5 

0.490 

O.OOJ 

U .490 

0 .2 68 0 

0.2866 

0.2862 

56/36. 

0.3334 

0. 3339 

6 

0.24C 

0 6 1 

0.335 

0.242C 

0.2255 

0.2436 

4Od50 . 

0.27C6 

0.2668 

7 

0. 300 


0. 3 86 

0.2550 

0.2350 

0.2591 

51359. 

0.2E71 

0.2956 

fc 

0. 360 


0.437 

0.25eC 

0.2353 

0 . 2 T 1 4 

55961. 

0.3C13 

0.3171 

9 

0.415 

0. 1 

0.483 

0.2860 

0.2299 

0.2830 

57122. 

C.3388 

0 .3342 

10 

0.2 40 

0 *5 J 

\J .404 

0.1850 

0. 2jJ0 

0. 1 796 

45613. 

U.2C74 

0.1994 

ll 

0. 3 JO 

06H^U 

O.503 

0.2050 

0.2072 

0.1803 

45253. 

0.2289 

0.2C78 

12 

C . 360 

G 6 5 c. 

J . 54 1 

0. 1640 

0. ^044 

0. 1 750 

4 84 8 7 . 

U .2064 

0 .2076 

13 

0.4 15 

0 •‘r 

0 .5 /6 

0 . 1820 

0. 1949 

0.1678 

48571. 

0.2237 

0.2 c30 

14 

0.2h0 


0.534 

0. 1500 

0. 1674 

0.1561 

35914. 

0.1770 

0. U32 

15 

0.3 00 

U • 

0. 5bo 

0. Io50 

0. 1 751 

0 . 1 6 j6 

3 7330. 

u . 1 6 70 

0 . 1 90 1 

\b 

'j . 3 6 0 

: .5^0 

J.598 

0.1330 

<^.1747 

0 . 1 5v; 7 

4.*208. 

0. 1 795 

0. 1664 

1 7 

0.4 15 

0. J 

0. 62 7 

0 . 1 45 C 

0.1685 

0.1418 

4 3 7 J 7 . 

U.1979 

0.1814 

ly 

0. 550 

0 

O.o98 

0. 1100 

0. 1368 

0.1157 

9ol9. 

J . 1 364 

0.1515 

10 

C . 1 80 

v.;20 

0.573 

0.125O 

0.1161 

0.1292 

2 7066. 

0.1642 

0. 1564 

?c 

0. 240 

0. 78^0 

0.599 

C. 1 iu u 

u. 1307 

0. I4l4 

31828. 

U . lo6 1 

0 • i 7 o3 

21 

0.3 00 

0 . 74.0 

0 . 6 2 4 

0. 1 300 

0. 1383 

0. 1454 

35534. 

0.1643 

0.1 864 

22 

0. 3oO 

0. 7: 0 

J . o 5 0 

0.1350 

J. 1400 

0.1432 

3ol38. 

u.l 7S7 

0.1883 

23 

0.4 15 

D. 74^0 

3.6 73 

0.1350 

o. 1372 

0.13 72 

3 71 70. 

U.I881 

0 . 1843 

24 

0. 1 95 

0 0 

J .6o5 

v^. 1280 

0.0771 

0.1069 

2U802. 

0.1732 

0.1378 

25 

0.2 85 

0# 6o0 

0.692 

0.132 0 

0.0886 

0.1226 

24640. 

0.1692 

0 . 1628 

2b 

0.3 75 

0# b80 

0. 71v 

0. 1 300 

0. 092 7 

0.1271 

25844. 

0.1655 

0 .1 730 

27 

C.475 

C.ddj 

0.749 

• 1 240 

0.0901 

0. 1222 

24C99, 

C.l 7C5 

0. 17ul 

28 

0.545 

0« ooO 

0. 769 

0.1 130 

0.0647 

0.1141 

24551. 

0.1577 

0.1608 

29 

0.6 20 

0 • b b 0 

J. 792 

0. 1120 

0. 07a5 

0.1026 

219ol . 

u;i617 

0.1456 
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in KSI 

B— 1.~ And nuclnun Allowabls buckling 

7075-T6 bar* *h**t and plat*. 











APPENDIX C 


STIFFNESS FORNULATIOM OP SUBSONIC 
STEADY-STATE NINO LOADING 

Thm purpos* of this sppsndix Is to sxplsln ths dsrlvstlon of ststlc 
squlllbrium conditions rsquirsd for ths calculation of ths subsonic stsady- 
stats wing loadings on a flsxibls wing. Ths dsrivation of ths asrodynamic 
principlss undsrlying this sppsndix ars suanarizsd hsrs and ars sxplainsd 
in dspth in rsfsrsncs 11. Ths assuaptions coonon to airfoil thsory that 
apply includs: (1) ths flow is potsntial with nsgligibls offsets from 

boundary laysrs, ssparation, or shocks; (2) ths wing thicknsss is sstall; 

(3) a stagnation point sxists at ths wing trailing sdgs; (4) ths angls of 
attack, a, is small so that tan a *• sin a ■ a (radians) ; and (5) ths span- 
wiss lift distribution dus to a givsn horssshos vortsx is constant ovsr its 
span of application. The aquations ussd to dsscribs ths stsady-stats squi- 
librium ars outlinsd bslow. 

Ths lift or circulation distribution, which variss along ths span of 
a wing, can bs visualizsd as resulting from a systsrn of horssshos vorticss, 
each of which is of constant strength. Such a system of horssshos vortices 
is illustrated in figures C-1 and C-2. It is obvious from ths figures that 
ths shape of ths actual load distribution may bs approximated to any degree 
of accuracy by a suitable number of horseshoe vortices. The relationship 
used to describe the lift of each vortex due to the local angle of attack 
at a given spanwise station is 

where [ is the aerodynamic- induction or dovmwash matrix, which contains 

the effects of wing planform geometry; q is the free-stream dynamic pres- 
sure; * matrix of t%«o-dimensional section lift-curve slopes; It) 

is a vector of section lifts per unit .span; and | is a vector of final 

section angles of attack. 

The final angle-of-attack variation across the span, |a^|,.'an be 
considered to be composed of four parts: 

l“fl * l“rl " l“ol * l“.l * l“«l 
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H«r« (<>f| th« angl« of Attack of th« root-soction soro-lift lino; 

I I la tho aoctlon anglo of attack duo to twiat | | and duo to doforma* 

tion roaultlng fro« tho proaonco of oxtornal atoroa or forcoa on tho wing 
i^cml * i^a} anglo of attack cauaod by atruetural dofloction of a 

floxiblo vlng duo to aorodynaaiic forcoa; and |a^| la tho Inducod anglo of 
attack duo to alloron or flap dofloction, A. 

Tho doformatlon porpondlcular to tho piano of tho wing la rolatod to 
tho aoctlon lift por unit apan aa followa: 

C’']K * ’'em) * 

whoro [k] la tho atlffnoaa matrix, *nd |w^| aro tho atructiiral 

doformatlona duo to aorodynamlc and concontratod maaa loadlnga roapoctlvoly, 
and lo tho tr ana formation rolatlng tho apanwlao loading to %fork 

oqulvalont loadlnga in tho atructuro. Thla tranafonaatlon la dotallod 
lator In thla appondlx. 

Tho romalnlng oquatlona doacrlblng tho atoady-otato oqulllbrlum con- 
dition aro dorlvod from tho atatlc balanco roqulromonta of tho vanlahlng 
of tho aummatlon of forcoa normal to tho wing and for tho vanlahlng of tho 
roaultlng momonta In pitch and roll. Tho forcoa and momonta roaultlng 
from tho aircraft body contor-of-gravlty load, W, tho tall load, P^, and 
tho oxtornally locatod concontratod maaaoa, ^ca* includod aa ahoim in 

figuroa C-1 and C-2. Tho offocta of roll aro includod in thla dorlvatlon 
by tho addition of a loading proportional to a control-aurfaco dofloction, 

< 5 . Whon 7ohiclo symmotry la aaaumod, all torma aaaoclatod with thla do- 
floction vanlah. Tho following aymmotry conatant ia dofinod: 

1 1, nonaymmotrlc planforma 
2, aymmetric planforma 

Summation of forcoa normal to tho wing glvoa 

e{2h|'^{f) - -P^ + nw + 0na»^ (C-4) 

Summation of momonta about tho pitch (y-) axla givoa 
e|2hx^)'^l/| - -P^^ + nwx^ ^ ^"(^em’^cm)* 
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Suonatlon of mononta about tha roll (x-) axla givao 


- na*^^ (c-6) 

Bara h la tha aaail-width of a glvan bound vortax, n la tha load factor, 
Xyt/v daflna tha location of tha cantrold of application of tha aactlon 
lift, x^,y^ daflna tha location of tha concantratad maaaaa, c la tha 


aactlon avaraga chord langth, and la tha rata of changa of aactlon- 

nomant coafflclant with raapact to control-aurfaca daflactlon. Whan ^ 1, 

panal coafflanta and maaaaa for both wlnga muat ba glvan. Zf 0 ■ 2, only 
tha aymmatrlc panal or maaa la uaad. 

Tha following aquatlona ralata tha componanta of aquation (C-2) to 
thalr approprlata varlablaa: 
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whara | >" 5 1 tha ratio of tha local Ilf t-curva slopa with raapact to 

control-aurfaca daflactlon, to tha local Ilf t-curva slopa with raapact to 
angla of attack and [*^( 2 ] !• tha tr ana formation ralatlng tha local aactlon 

angla of a. tack dua to daformatlon of tha corraapondlng structural dlsplaca- 
mants. Th:s transformation haa baan darlvad from a laast-squares fit of 
tha structural daformatlon In tha chordwisa direction at each section and 
is detailed later In this appendix. 

By use of tha aquatlona of thij appendix, the equilibrium conditions. 
Including pitch and roll trim, are derived as follows. From aquations (C-1) , 
(C-2), (C-3) , (C- 8 ) , and (C-9) , tha structural daformatlon Is expressed in 

terms of the root angle of attack, local twist, and control- surf ace deflec- 
tion : 

([k] - [t„]['4<I»J[s - w^j - » )a^j * ja,}) 


(C-11) 
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or 


1“. * ’'d • [’'k](l“t| * l“t| * hst) 
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Th« local section lift a^y than ba %nrittan as 


(‘) ^ I'^l * l“M * [''a]r'k]()“tl *• Kl * l“sl)] 


(C-13) 


From aquation (C-6) and tha results of (C-13), the control- surface deflec- 
tion required for trim in roll becooMs: 

For a symmetric planform (^ ■ 2), A ■ 0 

For a nonsymmetric planform O ■ 1) , 
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The root angle of attac)c is derived from equations (C-4) and C-S) and the 
results of (C-13) and (C-14) as 
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Th« tail load raqulrad for static aquilibriua bacomas, froa (C-4) , 


- nw + - pj2hj'^(/> (C-16) 

whara a positiva would rasult in a pitch-down momant for a raar- 

mountad tail configuration, and tha lift distribution, (I), is obta^nad 
by substitution of aquations (C-14) and (C-15) into (C-13). 


Intagration into WADES Program 

'Hia abryva systam of aquations has baan intagratad into tha WADES 
program to provide it with a simplifiad subsonic aarodynamic loads capa- 
bility. version is optionally available at load time during program 

axacution. Sas tha program documentation for usage 

A number of tha aarodynamic section properties are defaulted to 
theoretical vc«lues in the January 1975 version of tha program. Tha inter- 
nally defined properties are: (1) tha lift-curve slope * 6.28; 

(2) tha section-moment coefficient " 0; (3) no control- stir face 

momant, ]m^| ■ 0; and (4) tha aircraft body canter of gravity is sat 
approximately at tha aarodynamic canter. The loads analysis is restricted 
to symmetric configurations because of other program restrictions. To 
obtain a reasonable estimate of tha shape of tha spanwise load distribu- 
tion, the number of horseshoe vortices is internally limited to the number 
of degrees of freedom used to describe the normal deformation of tha wing 
surface. 
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O«riv«tion of Trans fonMtions and 

Tha transforvation [T^] as ussd in aquation (C-6) dafinss tha rala- 
tionship of tha structural daforvation to tha chanqa in local section angle 
of attack. Since tha WA0B8 prograai uses assumed sodas of deformation that 
allow for local chordwisa curvature, a relationship was derived to approx* 
imata tha rigid-body motion of tha section airfoil. Zn order to make this 
approximation, a laast-squaras procedure was derived to fit a linear func- 
tion to tlie deformed shape. The following equations outline that approxi- 
mation and the calculation of the resulting transformation. 

The least-squares fit is obtained by assuming that at the k‘th control 
point on the wing the structural deformation can be approximated by the 
linear function 


The least-squares problem then is defined as 

V - ^ 




(C-17) 


(C-18) 


where n is the nusdaer of chordwise points at which the function is susined. 
This problem may than be transformed to an integral form by multiplying 
both sides by the distance between points and taking the limit as the in- 
crement vanishes. Then (C-18) becomes 
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(C-19) 
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whara tha funetlona ^ daflnad i** taxms of tha axidal 

ahapaa of tha diaplacaaiant function aa 
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(C-21) 


Slnca tha angla of attack la defined aa tha chordwiaa alopa of tha function, 
tha angla dua to atructural dafornation can ba axpraaaad aa 

a 

^ Rde R 


(C-22) 


whara tha angla ia in radiana. By aolving for tha coafficiant a in 

9 

aquation (C-20) tha angla dua to atructural daformation at tha k'th control 
atation bacoaiaa 
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(C-23) 


Tha raaulting aat of aquationa uaing tha abova approximation avaluatad at 
aach of tha control pointa producaa tha ralationahip 

{a} • (C-24) 

Tha tranaformation dafinad in aquation (C-3) ralataa tha 

aaction lift to tha %#ork aquivalant loado ganaratad by applying tha lift 
aa a concantratad forca at ita appropriata control point. Tha work aquiv- 
alant load in tarma of tha variation of axcarnal %fork may ba writtan aa 
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•nd th«r«for« th« co«fflcl«nts of tho tranaformatlon may b« oxproasod a* 
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whara nondlaanalonal locatic'x* of tha control pointa 
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